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ABSTRACT 

We present an end-to-end, two-phase model for the origin of globular clusters (GCs). In 
the model, populations of stellar clusters form in the high-pressure discs of high-redshift 
(z > 2) galaxies (a rapid-disruption phase due to tidal perturbations from the dense interstel¬ 
lar medium), after which the galaxy mergers associated with hierarchical galaxy formation 
redistribute the surviving, massive clusters into the galaxy haloes, where they remain until 
the present day (a slow-disruption phase due to tidal evaporation). The high galaxy merger 
rates of ^ > 2 galaxies allow these clusters to be ‘liberated’ into the galaxy haloes before 
they are disrupted within the high-density discs. This physically-motivated toy model is the 
first to include the rapid-disruption phase, which is shown to be essential for simultaneously 
reproducing the wide variety of properties of observed GC systems, such as their universal 
characteristic mass-scale, the dependence of the specific frequency on metallicity and galaxy 
mass, the GC system mass-halo mass relation, the constant number of GCs per unit supermas- 
sive black hole mass, and the colour bimodality of GC systems. The model predicts that most 
of these observables were already in place at z = 1-2, although under rare circumstances GCs 
may still form in present-day galaxies. In addition, the model provides important constraints 
on models for multiple stellar populations in GCs by putting limits on initial GC masses and 
the amount of pristine gas accretion. The paper is concluded with a discussion of these and 
several other predictions and implications, as well as the main open questions in the field. 

Key words: galaxies: evolution — galaxies: formation — galaxies: haloes — galaxies: star 
clusters - globular clusters: general — stars: formation 


1 INTRODUCTION 

All massive galaxies (Mstar > 10^ M©) in the local Universe 
host populations of globular clusters (GCs), gravitationally bound 
stellar systems that are typically old (r ~ 10^° yr) and massive 
(M ~ 10^ M©). The origin of GCs is a major unsolved problem 
on the interface between star and galaxy formation. In part, our un¬ 
derstanding of GC formation is limited because most GCs in the 
Universe must have formed at redshifts z > 2. At such distances, 
the sizes of star-forming giant molecular clouds (GMCs) are gen¬ 
erally unresolved, implying that the physical processes leading to 
GC formation must be inferred indirectly. Important, independent 
examples of such constraints on the GC formation process exist 
from the formation of massive clusters in the local Universe, the 
observed conditions for star and cluster formation in high-redshift 
galaxies, and the present-day properties of GC populations. 

In this work, we present a new model for GC formation that 
combines the current observational and theoretical constraints. This 
model is the first to provide an end-to-end (yet simple) descrip¬ 
tion for the origin of GCs, from their formation at high redshift, 
through the relevant evolutionary processes, until the present day. 
It explains the observed properties of GC populations as the natural 


outcome of regular star and cluster formation in the high-redshift 
Universe. 

Before discussing the origin of GCs further, we should de¬ 
fine the term ‘globular cluster’. In the literature, a large variety 
of definitions has been used, based on metallicity (‘metal-poor’), 
mass (M = 10^-10® M©), age (r ~ 10^° yr), location (‘in the 
halo’), or chemical abundance patterns (‘multiple stellar popula¬ 
tions’). However, exceptions to all of these criteria can be given. 
For instance, GCs are known to exist with (super-)solar metallici- 
ties, with masses of M < 10^ M©, with ages r ~ 5 Gyr (or even 
as young as a few 100 Myr, Schweizer & Seitzer 1998), with kine¬ 
matics and positions consistent with the host galaxy’s bulge com¬ 
ponent, or without multiple stellar populations (e.g. Harris 1996; 
Dinescu, Girard & van Altena 1999; Carraro et al. 2006; Forbes & 
Bridges 2010; Walker et al. 2011). In this work, we will therefore 
use a definition which will be shown to be physically more fun¬ 
damental. The term ‘globular cluster’ will be used to refer to any 
object satisfying the following condition. 

A gravitationally-bound, stellar cluster that in terms of its po¬ 
sition and velocity vectors does not coincide with the presently star¬ 
forming component of its host galaxy. 

In this definition, the term ‘stellar cluster’ excludes dark 
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matter-dominated systems such as dwarf galaxies, and ‘the 
presently star-forming component’ refers to the part of position- 
velocity space occupied by a galaxy’s star-forming gas. Of course, 
a GC may pass through this component - in this case, it coincides 
in position, but not in velocity space. It may also co-rotate with the 
star-forming component, but in a position above or below the galac¬ 
tic plane (as holds for part of the GC population of the Milky Way, 
see Frenk & White 1980). We have deliberately omitted limits on 
the GC metallicity, mass, or age. However, it is one of the goals of 
this paper to show that the ranges of these quantities observed in 
GCs are implied by a combination of the above definition and the 
physics of cluster formation and evolution. 

Early theories aiming to explain the origin of GCs often in¬ 
voked the special conditions in the early Universe (e.g. high Jeans 
masses prior to the formation of the first galaxies or in young galaxy 
haloes, see Peebles & Dicke 1968; Fall & Rees 1985). More re¬ 
cent models have begun to focus on the dynamical evolution of 
GCs within their host galaxy haloes in an attempt to reconstruct 
their present-day properties from an assumed initial population of 
stellar clusters (e.g. Fall & Zhang 2001; Prieto & Gnedin 2008; 
Tonini 2013; Li & Gnedin 2014). These theories generally omit a 
description of the actual GC formation process itself and may there¬ 
fore be missing important physics. Another family of models aims 
to constrain GC formation physics using the chemical abundance 
variations observed in present-day GCs (e.g. Decressin et al. 2007; 
D’Ercole et al. 2008; Conroy & Spergel 2011; Krause et al. 2013; 
Bastian et al. 2013b). These models often require new (and thus far 
unobserved, see Bastian et al. 2013a; Kruijssen 2014) physics to 
explain the observed abundance variations. 

In most of the above cases, the formation of GCs is reverse- 
engineered using their present-day properties and/or their recent 
dynamical evolution. However, further constraints can be added by 
considering the formation environment of GCs. The old ages of 
most GCs imply that the majority of them must have formed at 
z > 2, close to the peak of the cosmic star formation history at 
z = 2-3 (e.g. Hopkins & Beacom 2006). Indeed, observations of 
galaxies at these redshifts reveal vigorously star-forming systems, 
with star formation rates (SFRs) exceeding several 10^ M© yr“^ 
(e.g. Forster Schreiber et al. 2009; Daddi et al. 2010). The gas con¬ 
tent of these high-redshift galaxies has been surveyed in great de¬ 
tail over the past couple of years (e.g. Genzel et al. 2010; Tacconi 
et al. 2010, 2013; Swinbank et al. 2011), revealing high gas densi¬ 
ties (E = 10^-10^ Mq pc“^) and turbulent velocities (a = 10- 
100 km s“^), implying gas pressures (P/k ~ 10^ K cm“^) some 
three orders of magnitude higher than in the discs of nearby galax¬ 
ies such as the Milky Way (e.g. Kruijssen & Longmore 2013). This 
typical pressure falls nicely in the range required for GC formation 
(P/k = 10^-10® K cm“^, see Elmegreen & Efremov 1997). 

While the high-pressure conditions of star formation in z > 2 
galaxies may seem unique, analogous environments exist in the lo¬ 
cal Universe. Galaxy mergers, starburst dwarf galaxies, and galactic 
nuclei are all known to reach the extreme pressures seen regularly 
at high redshift (e.g. Downes & Solomon 1998; Kruijssen & Long- 
more 2013). Such high-pressure environments are almost univer¬ 
sally seen to produce stellar clusters that are indistinguishable from 
GCs in terms of their masses and radii (Portegies Zwart, McMil¬ 
lan & Gieles 2010; Longmore et al. 2014; Kruijssen 2014). These 
‘young massive clusters’ (YMCs; with masses M = 10^-10® M©, 
radii rn = 0.5-10 pc, and ages r < 1 Gyr) had already been 
discovered in the 1980s (Schweizer 1982, 1987), but their number 
increased spectacularly after the launch of the Hubble Space Tele¬ 
scope (HST; see e.g. Holtzman et al. 1992; Schweizer et al. 1996; 


Whitmore et al. 1999; Bastian et al. 2006). These observations un¬ 
equivocally showed that GC-like clusters are still forming in the 
local Universe wherever the conditions are similar to those seen in 
star-forming galaxies at the peak of the cosmic star formation his¬ 
tory at z = 2-3. The obvious difference between YMCs and GCs 
is that the former are generally still associated with the presently 
star-forming components of their host galaxies (thereby disquali¬ 
fying them as GCs in the above definition), whereas the latter are 
not. In this paper, it is proposed that this separation occurs naturally 
during galaxy mergers and the accretion of satellite galaxies, both 
of which are common events in the context of hierarchical galaxy 
formation (White & Rees 1978; White & Frenk 1991). 

Given the ubiquity of GC-like cluster formation under high- 
pressure conditions, it is undesirable to devise ‘special’ physical 
mechanisms to explain the origin of GCs. Instead, the primary 
question(s) to ask should be: 

(i) Could the products of regular cluster formation in high- 
redshift galaxies have survived until the present day? 

(ii) If so, are these relics consistent with the properties of 
present-day GC populations? 

It is the aim of this work to address the above questions by con¬ 
structing a simple model for (1) the formation of regular stellar 
clusters in the early Universe and (2) their survival until the present- 
day in the context of hierarchical galaxy formation. To do so, we 
combine the current observational and theoretical constraints on 
nearby YMC formation, high-redshift star formation, and present- 
day GC populations. While no fundamentally new physical mech¬ 
anisms are added, this is the first time that an end-to-end model 
for the cosmological formation and evolution of massive clusters 
includes a physical description for the formation and early evolu¬ 
tion of these clusters within the high-pressure environments of their 
host galaxies. We show that this initial phase plays an essential role 
in setting the properties of the surviving cluster population (its im¬ 
portance was also highlighted by Elmegreen 2010 and Kruijssen 
et al. 2012b). As will be discussed at length, the modelled stellar 
cluster populations at z = 0 match a wide range of independent, 
observed properties of present-day GC populations, such as the GC 
mass function, specific frequency, and the constant GC population 
mass per unit dark matter halo mass. 

Of course, this is not the first paper aiming to put GCs in the 
context of galaxy formation. However, previous efforts to explain 
the properties of GC populations have focused on: 

(i) exploring special, high-redshift conditions in which GCs 
could conceivably have formed (e.g. in protogalactic clouds with 
high Jeans masses, in galaxy mergers, or during reionisation, see 
Peebles & Dicke 1968; Fall & Rees 1985; Ashman & Zepf 1992; 
Katz & Ricotti 2014); 

(ii) modelling a Hubble time of dynamical evolution and evap¬ 
oration of GCs in the haloes where they presently reside, with the 
aim of explaining the differences between observed young cluster 
populations and old GCs (e.g. Prieto & Gnedin 2008; Muratov & 
Gnedin 2010; Katz & Ricotti 2014; Li & Gnedin 2014); 

(iii) following the GC population during hierarchical galaxy for¬ 
mation by assuming that some fraction of high-redshift star forma¬ 
tion occurred in the form of GCs, without accounting for any cluster 
disruption (e.g. Kravtsov & Gnedin 2005; Tonini 2013). 

The theory in this paper differs fundamentally from these previ¬ 
ous approaches. Instead of exclusively focusing on special condi¬ 
tions for GC formation, disruption by tidal evaporation, or hierar¬ 
chical galaxy growth, in this paper we take the “Occam’s Razor” 
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approach of combining all mechanisms that are known to affect the 
formation and evolution of stellar clusters, and to put these in the 
context of galaxy formation to verify if the properties of observed 
GC populations are retrieved. The resulting model is the first to 
include the environmental dependence of GC formation and early 
disruption by tidal shocks in the host galaxy disc. As mentioned 
above, we will show that this early phase is dominant in setting 
the properties of the z = 0 GC population. In addition, we will 
improve on previous models for GC evaporation in galaxy haloes 
(which used the environmentally-independent expression of Spitzer 
1987) by following recent A^-body simulations (e.g. Baumgardt & 
Makino 2003; Gieles & Baumgardt 2008; Lamers, Baumgardt & 
Gieles 2010) and accounting for the environmental dependence of 
GC evaporation due to the presence of a tidal field. 

The paper is structured as follows. In §2, we will first sum¬ 
marise the main observational (§2.1) and theoretical (§2.2) con¬ 
straints on GC formation. In §3, we will combine these constraints 
in a simple, two-phase model for GC formation in which the vio¬ 
lent conditions at high redshift promote GC formation due to (1) 
the high initial cluster masses that are attainable in high-pressure 
environments and (2) the rapid migration of (Y)MCs into the qui¬ 
escent galaxy halo due to the high incidence of galaxy mergers at 
z > 2, both of which enable the long-term survival of GCs. The 
predictions of the model are contrasted with a wide range of ob¬ 
servational data, illustrating how the observed properties of the GC 
population have been sculpted by the violent conditions of their 
formation (§4). We also discuss the implications for our current 
understanding of GC formation and provide predictions for future 
work (§5). The main open questions in the field are discussed in §6, 
and the summary and conclusions of this work are provided in §7. 


2 CURRENT CONSTRAINTS ON GLOBULAR CLUSTER 
EORMATION 

In this section, we summarise the constraints on GC formation (also 
see the recent review of Kruijssen 2014). These constraints are 
drawn from the fields of nearby YMC formation, high-redshift star 
formation, and present-day GC populations, as well as the theory of 
cluster formation and evolution. Our new model will be described 
in §3 and is based on the constraints discussed in this summary. 
The comparison to observations in §4 also focuses on the observ¬ 
ables described here. The reader who is exclusively interested in 
the model is advised to skip this section and continue to §3 and §4. 

2.1 Observational constraints 

2.1.1 YMCs in the local Universe as analogues for GC formation 

Ever since the discovery of YMCs with masses similar to (or larger 
than) those of GCs in the 1980s and 1990s, numerous authors have 
suggested that these YMCs could be the progenitors of future GCs. 
Initially, the formation of these extreme stellar clusters was thought 
to be exclusive to galaxy mergers or merger remnants (e.g. Ash¬ 
man & Zepf 1992), but YMCs have now been seen in a variety 
of high-pressure environments, including starburst dwarf galaxies 
and galaxy centres (e.g. Eiger, McLean & Morris 1999; Eiger et al. 
2002; Anders et al. 2004). If GCs indeed formed as regular YMCs, 
then the characterisation of the YMC formation process in the local 
Universe provides insight into the formation of GCs. 

There is a wide range of literature on YMC formation (see e.g. 
Portegies Zwart, McMillan & Gieles 2010; Longmore et al. 2014). 


The observational constraints most relevant in the context of GC 
formation are summarised here. 

(i) YMCs are seen to form through the hierarchical merging of 
smaller structures (Elmegreen & Elmegreen 2001; Bastian et al. 
2007; Sabbi et al. 2012; Gouliermis, Hony & Klessen 2014; Rath- 
borne et al. 2015; Walker et al. 2015) on a short (~ Myr) time-scale 
(Elmegreen 2000; Bastian et al. 2013a; Cabrera-Ziri et al. 2014). 
After this time, no gas is left within the YMC due to a combination 
of gas consumption in star formation and gas expulsion by stel¬ 
lar feedback and the YMC resides in virial equilibrium (Rochau 
et al. 2010; Cottaar et al. 2012; Henault-Brunet et al. 2012; Clark¬ 
son et al. 2012). 

(ii) Only some fraction of all star formation occurs in bound 
stellar clusters (Lada & Lada 2003), which is referred to as the the 
cluster formation efficiency (CEE or L, Bastian 2008). Within the 
hierarchical density structure of the interstellar medium (ISM), the 
highest-density peaks achieve the highest star formation efficien¬ 
cies (SFEs or e), allowing them to remain gravitationally bound 
when star formation is halted by feedback (Elmegreen & Efremov 
1997; Elmegreen 2008; Kruijssen 2012; Wright et al. 2014). Be¬ 
cause the mid-plane gas density increases with the gas pressure, 
the CEE is seen to increase with the pressure too, from L ~ 0.01 at 
P/k < 10^ K cm“^ to r ~ 0.5 at P/k > 10^ K cm“^ (Larsen 
2000; Goddard, Bastian & Kennicutt 2010; Cook et al. 2012; Silva- 
Villa, Adamo & Bastian 2013; Adamo et al. 2015). 

(iii) The initial cluster mass function (ICMF) follows a power 
law (diV/dMi oc Mf) with index a = —2 for cluster masses M > 
10^ M© (Zhang & Fall 1999; Lada & Lada 2003; Larsen 2009). At 
the high-mass end, the ICMF has a truncation that is well-fit by 
an exponential function [dN/dMi (x Mf exp (—Mi/Me), where 
Me is the truncation mass, see e.g. Gieles et al. 2006a; Larsen 2009; 
Bastian et al. 2012], which was originally proposed by Schechter 
(1976) to describe the galaxy mass function. 

(iv) The ICMF truncation mass is proportional to the Toomre 

(1964) mass Mt = of the host galaxy disc (Kruijssen 

2014; Adamo et al. 2015), where a is the gas velocity dispersion 
and S is the gas surface density. This mass-scale represents the two- 
dimensional Jeans mass, i.e. the maximum mass-scale below which 
self-gravity can overcome galactic shear (Elmegreen 1983; Kim & 
Ostriker 2001), and determines the mass of the most massive GMC 
in a galaxy. The mass of the most massive cluster Mt,c 1 is subse¬ 
quently obtained by multiplying the Toomre mass by the SEE and 
the CEE, i.e. Mt,c 1 = cTMt. The Toomre mass steeply increases 
with the pressure, hence the mass of the most massive cluster in a 
galaxy is a good proxy for the pressure in its formation environ¬ 
ment. Indeed, nearby galaxies that produce clusters with masses 
similar to the most massive GCs (M >10® M©) have ISM pres¬ 
sures similar to those seen at high redshift (Kruijssen 2014). 

(v) YMCs have radii (rn = 0.5-10 pc) that are roughly inde¬ 
pendent of their masses (Larsen 2004). The variation of the initial 
cluster radius with the galactic environment is an open problem - 
a slight decrease of the radius with the ambient density may ex¬ 
ist. This would be expected if clusters born in strong tidal fields 
(or high-pressure environments) are more compact than those in 
weak tidal fields (e.g. Elmegreen 2008). Scaling the typical YMC 
radius in the Galactic disc to the strong tidal field of the Galactic 
centre does reproduce the small (0.5 pc) radius of the Arches clus¬ 
ter, which resides within the central ~ 100 pc of the Milky Way. 
While this comparison is certainly suggestive, this is an area where 
current constraints must be improved. 

The above points show that the formation of GC-like YMCs (M ^ 
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10® Mq, see §4 and §5) occurs naturally under high-pressure con- 
ditions, due to the elevated CFEs and Toomre masses in such envi¬ 
ronments. 


2.7.2 Star formation in high-redshift galaxies as the 
environmental conditions for GC formation 

Most GCs have ages in the range r = 10-13 Gyr (e.g. Forbes & 
Bridges 2010; VandenBerg et al. 2013), indicating that the peak of 
cosmic GC formation must occur at redshift z = 2-6. In this paper, 
a fiducial GC formation redshift of z ~ 3 is adopted.^ Galaxies at 
these redshifts are too distant to observe the GC formation process 
in detail, but the global conditions for star formation in z > 2 
galaxies (i.e. near the peak of the cosmic star formation history at 
2 : = 2-3, Hopkins & Beacom 2006) can provide further constraints 
on the formation of GCs (e.g. Shapiro, Genzel & Forster Schreiber 
2010; Elmegreen 2010; Kruijssen et al. 2012b). For instance, they 
can help answer whether the properties of the ISM in these galaxies 
promote the formation of GC-like YMCs. 

Like for YMC formation, there exists a large body of literature 
on star formation in high-redshift galaxies. Here, we summarise the 
properties of z > 2 galaxies relevant for GC formation. 


(i) High-redshift galaxies with stellar masses Mstar > 4 x 
10® M© are characterised by high gas fractions (/gas ~ 0.5, 
Tacconi et al. 2013) and specific star formation rates (sSFR = 
SFR/Mstar ~ 3 Gyr“^, Bouche et al. 2010; Rodighiero et al. 
2011), indicating that a large fraction of their stellar mass is in the 
process of being formed. The star formation rate per unit gas mass 
(or alternatively the gas depletion time tdepi = SFR/Mgas) is sim¬ 
ilar to that in nearby galaxies (tdepi 10 ® yr, Bigiel et al. 2011; 
Schruba et al. 2011; Tacconi et al. 2013; Fisher et al. 2014). 

(ii) The metallicities of star-forming galaxies atz^S are lower 
than those in the local Universe by a factor of ~ 5, but the overall 
shape of the galaxy mass-metallicity relation is similar. In Kruijs¬ 
sen (2014), we parameterised the mass-metallicity relations of Erb 
et al. (2006) and Mannucci et al. (2009) as 


[Fe/H] 
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(iv) In particular, the gas-rich discs of z > 2 galaxies exhibit 
clumpy morphologies (Elmegreen & Elmegreen 2005; Genzel et al. 
2011; Tacconi et al. 2013), suggesting widespread gravitational in¬ 
stabilities. These clumps have size-scales of ~ 1 kpc and masses 
M ^ 10^ M©, and may provide a natural formation environment 
for massive stellar clusters (Shapiro, Genzel & Forster Schreiber 
2010). Note that most gravitational instabilities form at large galac- 
tocentric radii (out to several kpc, e.g. Swinbank et al. 2012), im¬ 
plying that even the most massive clusters forming within them 
would require > 10^® yr to spiral in to the galaxy centre by dy¬ 
namical friction. 

(v) While the high-pressure conditions in gas-rich galaxies pro¬ 
mote the formation of massive clusters, the correspondingly-high 
gas densities should also limit the long-term survival of stellar clus¬ 
ters due to the frequent and strong tidal perturbations by dense gas 
pockets. These ‘tidal shocks’ have been shown to be the dominant 
cluster disruption agent in gas-rich environments, both in empiri¬ 
cal studies (Earners, Gieles & Portegies Zwart 2005) and theoreti¬ 
cal work (e.g. Earners & Gieles 2006; Elmegreen & Hunter 2010; 
Kruijssen et al. 201 1). Extrapolating the empirically-estimated life¬ 
times of clusters in local-Universe, high-pressure conditions (e.g. 
Gieles et al. 2005) suggest that M ^ 10® M© clusters would sur¬ 
vive for tdis ^ 10® yr in z > 2 galaxies. 

We see that even though massive cluster formation should take 
place under the high-pressure conditions of z > 2 galaxies, the 
lifetimes of stellar clusters in such environments are very limited. 
In order to survive for a Hubble time, the GC progenitor clusters 
must therefore have escaped the gas-rich bodies of their host galax¬ 
ies before they got disrupted (see §2.2 below). 

2.1.3 Present-day GC populations as the outcome of the process 
we are seeking to understand 

The product of the GC formation process and a Hubble time of GC 
evolution is readily observable in nearby galaxies. Any theory of 
GC formation should yield observables at z = 0 consistent with 
the properties of present-day GC populations. There is a wealth of 
reviews discussing GCs in the nearby Universe (e.g. Harris 1991; 
Brodie & Strader 2006; Gratton, Carretta & Bragaglia 2012). The 
observables most relevant to this work are summarised below, again 
drawing from the most recent review of Kruijssen (2014). 


for the redshift range 2 ; = 2-4. This functional form suggests 
that at z ~ 3, clusters with typical GC metallicities (i.e. between 
[Fe/H] ~ —1.6 and [Fe/H] ~ —0.5, Peng et al. 2006) can form in 
galaxies with stellar masses Mstar = 10^-10^^ M©. 

(iii) In high-redshift galaxies, the formation of stars and stel¬ 
lar clusters proceeds at high densities and pressures, with gas sur¬ 
face densities S = 10^-10® ® M© pc“^, SFR surface densities 
SsFR = 10“^-10® ® M© yr ^ kpc and gas velocity dis¬ 
persions (7 = 10-100 kms“^ (Genzel et al. 2010; Swinbank 
et al. 2012; Tacconi et al. 2013), implying pressures of P/k > 
10^ K cm“®. These conditions give rise to high CFEs (F ~ 0.5, 
Kruijssen 2012) and Toomre masses (Mt ^10® M© and there¬ 
fore Mt,c 1 ^ 10® ®, Kruijssen 2014), both of which promote the 
formation of GC-like YMCs. 

^ We emphasise that the formation redshift of 2 ; ~ 3 is not rigid and we 
certainly do not mean to suggest that GC formation is restricted to 2 ; = 2-6. 
It is only intended to represent a reasonable choice for the majority of the 
GC population. In fact, one of the main conclusions of this work will be 
that GCs are still forming in the present-day Universe. 


(i) The globular cluster mass function (GCMF) is peaked, with 
a near-universal peak mass of Mpeak ~ 2 x 10® M© (Jordan et al. 
2007), in strong contrast with the Schechter-type ICMF seen for 
young stellar cluster populations (see §2.1.1). Only in galaxies with 
present-day stellar masses Mstar < 10^® M©, a slight decrease of 
the peak mass with decreasing galaxy mass is observed (Jordan 
et al. 2007). This weak trend is currently unexplained. 

(ii) The near-universal peak mass is often interpreted as the re¬ 
sult of GC evaporation, which leads to the gradual disruption of 
low-mass GCs over a Hubble time (Vesperini 2001; Fall & Zhang 
2001), thereby turning the Schechter-type ICMF into the peaked 
GCMF. However, in this scenario the peak mass should vary both 
with the host galaxy and the location within the halo, because the 
evaporation rate increases linearly with the orbital angular velocity 
(see §2.2 and e.g. Vesperini & Heggie 1997; Portegies Zwart et al. 
1998; Baumgardt & Makino 2003). Such a strong environmental 
variation is is not observed (Vesperini et al. 2003; Jordan et al. 
2007; McLaughlin & Fall 2008). This calls into question the im¬ 
portance of evaporation in setting the peaked shape of the present- 
day GCMF. Clearly, an end-to-end understanding of the emergence 
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of the GCMF from the ICMF cannot be achieved without more de¬ 
tailed knowledge of the mass loss histories of GCs. 

(iii) The GCMF has an exponential truncation at the high-mass 
end (Fall & Zhang 2001; Kruijssen & Portegies Zwart 2009), which 
depends strongly on the host galaxy mass. In Milky Way-mass 
galaxies the typical truncation mass is Me ~ 3 x 10® M©, but 
it is an order of magnitude lower in Mstar ~ 10^ M© galaxies 
(Jordan et al. 2007). Considering that the truncation mass of the 
ICMF of young stellar cluster populations is proportional to the 
Toomre mass (see §2.1.1, Kruijssen 2014 and Adamo et al. 2015), it 
is plausible that this variation reflects a trend of increasing Toomre 
mass with galaxy mass (and metallicity, see § 2 . 1 . 2 ) at the time of 
GC formation. 

(iv) The importance of GC evaporation for the present-day 
statistics of GC populations is called into question further by con¬ 
sidering the number of GCs per unit stellar mass in the host galaxy 
[i.e. the specifle frequency Tjv = Ngc/(10^ M©)] as a function 
of the galactocentric radius. It is shown in Figure 1 that the spe¬ 
cifle frequency strongly decreases with increasing metallicity (and 
hence the host galaxy mass at formation, see §2.1.2). This may 
result in extemely high specifle frequencies in metal-poor dwarf 
galaxies like the Fornax dSph, IKN, and WLM, where up to 20% 
of low-metallicity ([Fe/H] < — 2) stars can reside in GCs (Larsen, 
Strader & Brodie 2012; Larsen et al. 2014). However, at flxed 
metallicity, there is no dependence of the specifle frequency on the 
present-day galactocentric radius (Harris & Harris 2002; Lamers 
et al. 2015). This suggests a minor impact of evaporation (which 
would have required a clear radial variation of Tjv, see § 2 . 2 ) and 
challenges models aiming to explain the present-day properties of 
the GC population as the result of evaporation (e.g. Fall & Zhang 
2001; Li & Gnedin 2014). Instead, it seems more likely that the 
specifle frequency was set at the time of GC formation or during 
the early evolution of GCs, before they reached their present spa¬ 
tial conflguration within galaxy haloes. In that case, the relation be¬ 
tween Tjv and [Fe/H] (or the host galaxy mass at formation) may 
be fundamental in setting other observed correlations between the 
properties of GC systems and their host galaxies (see below). 

(v) The total mass of a galaxy’s GC population is a near¬ 
constant fraction of the galaxy’s dark matter halo mass, i.e. log 77 = 
log (MGc,tot/Mh) = —4.5 in the halo mass range Mh = 10®- 
10^® M© (e.g. Spider & Forbes 2009; Georgiev et al. 2010; Harris, 
Harris & Alessi 2013; Hudson, Harris & Harris 2014; Durrell et al. 
2014). It is unclear if this linear relation between GC system mass 
and halo mass reflects a fundamental connection between dark mat¬ 
ter and GCs. Even if such a relation initially existed, it must have 
been affected by a Hubble time of GC disruption. 

(vi) The GC colour and metallicity distributions are bimodal 
(e.g. Zinn 1985; Peng et al. 2006). While the spatial distribution and 
kinematics of the metal-rich GC sub-population bear some imprint 
of the host galaxy’s spheroid, the metal-poor GCs are mainly asso¬ 
ciated with the stellar halo. Recent theoretical work has shown that 
the metallicity bimodality of GCs may emerge naturally from hi¬ 
erarchical galaxy formation (Tonini 2013). In this scenario, metal- 
rich GCs are formed together with the main spheroid of the galaxy 
and metal-poor GCs accreted by the tidal stripping of satellite 
(dwarf) galaxies. However, this model uses the z = 0 relation 
between galaxy mass and specifle frequency (which is intimately 
related to the relation of Figure 1) to initialise the GC population 
at the time of GC formation. It therefore omits the physics of GC 
formation and disruption. 
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Figure 1. This figure shows that the specific frequency of GCs does not de¬ 
pend on the present galactocentric radius and therefore must have been set 
during their formation or early evolution. Shown is the specific frequency 
Tjv as a function of metallicity [Fe/H] for GCs in the inner region of 
NGC 5128 (triangles, 10 kpc; data from Harris & Harris 2002), the 
outer region of NGC 5128 (diamonds, R ^ 10 kpc; data from Harris & 
Harris 2002), and in the Fornax dSph (square; data from Larsen, Strader & 
Brodie 2012). The division into radial bins for the GCs in NGC 5128 shows 
that Tjv chiefiy depends on [Fe/H], exhibiting little variation with the ra¬ 
dius at fixed [Fe/H]. The data point for the Fornax GC system is included 
solely to illustrate the metallicity trend to lower metallicities. 


2.7.4 Chemical abundance patterns 

Next to the characteristics of the GC population, additional con¬ 
straints on GC formation may potentially be obtained from their 
composition. GCs are observed to host multiple stellar populations 
in terms of their light element abundances (e.g. Gratton, Carretta 
& Bragaglia 2012), such as an anti-correlation between the Na and 
O abundances of their constituent stars. This anti-correlation is al¬ 
most a defining feature of GCs, as it has been found in all consid¬ 
ered GCs (with masses M ^ 3 x 10^ M©, e.g. Gratton, Carretta 
& Bragaglia 2012). These abundance patterns have recently been 
heavily exploited to construct scenarios for the formation of GCs, 
almost exclusively from a stellar evolutionary perspective. 

Three main models aiming to explain the observed abundance 
patterns exist (see Kruijssen 2014 for a detailed discussion), but 
each still have substantial problems to overcome - no definitive 
model for the multiple stellar populations in GCs has been achieved 
(Bastian, Cabrera-Ziri & Salaris 2015). The viability of the exist¬ 
ing models should be tested by comparing their predictions and 
assumptions to independent constraints on GC formation, i.e. by 
not only considering stellar-evolutionary constraints, but also in¬ 
cluding what is known about GC formation and evolution. Here, 
we summarise the chief models for later reference in §5. 

(i) Enrichment by winds from asymptotic giant branch (AGB) 
stars (e.g. D’Ercole et al. 2008; Conroy & Spergel 2011). This 
model reproduces the abundance variations with the formation 
of a second generation of stars after enrichment by AGB ejecta 
(> 30 Myr after the formation of the first generation). 

(ii) Enrichment by winds from fast-rotating massive stars 
(ERMSs; e.g. Decressin et al. 2007; Krause et al. 2013). This model 
reproduces the abundance variations with the formation of a second 
generation of stars from retained ERMS ejecta (< 10 Myr after the 
formation of the first generation). 
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(iii) Enrichment by the early disc accretion (EDA) of winds 
from FRMSs and massive interacting binaries (MIBs; de Mink et al. 
2009; Bastian et al. 2013b). This model reproduces the abundance 
variations by the sweeping up of enriched ejecta from ERMSs and 
MIBs with the protoplanetary discs around pre-existing low-mass 
stars. Out of the three models discussed here, it is the only model 
that does not require the formation of a second generation of stars. 

The AGB and ERMS models require GCs to have been 10-100 
times more massive at birth to create an enriched gas reservoir 
sufficiently massive reproduce to the observed, 1:1 ratio between 
enriched and unenriched stars (the mass budget problem). In addi¬ 
tion, the AGB model requires substantial pristine gas accretion to 
turn the abundance patterns resulting from stellar evolution into the 
observed ones (the gas accretion problem). By contrast, the ERMS 
model needs to retain gas for ~ 30 Myr after the formation of the 
first generation despite vigorous feedback from massive stars (the 
gas retention problem). Finally, the EDA model requires protoplan¬ 
etary discs to survive for 5-10 Myr (the disc lifetime problem), 
which may be longer than discs survive in nearby star-forming re¬ 
gions (Haisch, Lada & Lada 2001; Bell et al. 2013), especially at 
the extreme stellar densities of GCs (Adams et al. 2004; de Juan 
Ovelar et al. 2012; Rosotti et al. 2014). 

We will revisit these problems in §5, by addressing questions 
such as how much more massive GCs could have been at birth given 
their disruption histories, or how much mass young GCs could have 
accreted under the conditions of high-redshift galaxies. 

2.2 Physics of globular cluster formation and evolution 

We now briefiy summarise the physical processes governing the 
formation of GCs and their evolution until the present day. 

(i) Whereas this work has the specific aim of testing whether 
GCs may be the descendants of regular YMC formation at high 
redshift, other GC formation theories exist too. It was proposed that 
GC from in mergers (Ashman & Zepf 1992), but it is now clear that 
the universally-high Toomre mass in the z > 2 Universe enables 
the formation of GC-like YMCs in normal disc galaxies too (e.g. 
Shapiro, Genzel & Forster Schreiber 2010). Likewise, recent work 
(e.g. Forbes & Bridges 2010; Georgiev et al. 2010; Conroy, Loeb 
& Spergel 2011) seems to rule out that GC formation is the result 
of a high Jeans mass following recombination (Peebles & Dicke 
1968), thermal instabilities in hot galaxy haloes (Fall & Rees 1985), 
or star formation before reionization within individual dark matter 
haloes (e.g. Peebles 1984; Bekki 2006). However, a small fraction 
of the GC population may represent the former nuclei of tidally 
stripped dwarf galaxies (Mackey & van den Bergh 2005; Lee et al. 
2009; Hartmann et al. 2011). The fraction of GCs formed this way 
is likely small, i.e. less than 15% (Kruijssen & Cooper 2012). 

(ii) Numerical simulations of YMC formation have advanced 
greatly in the last decade (see Kruijssen 2013 for a recent review). 
These models now reproduce the observational picture of the hi¬ 
erarchical fragmentation of the ISM (e.g. Tilley & Pudritz 2007; 
Bonnell, Clark & Bate 2008; Krumholz, Klein & McKee 2012) 
and the merging of stellar aggregates (Maschberger et al. 2010; Fu¬ 
ji!, Saitoh & Portegies Zwart 2012), in which gravitationally bound 
structure naturally arises at the highest-density peaks (e.g. Offner 
et al. 2009; Kruijssen et al. 2012a; Girichidis et al. 2012). On a 
galactic scale, this model of cluster formation in high-density peaks 
within a hierarchical ISM naturally reproduces the observation that 
the fraction of star formation occurring in bound clusters (the CEE) 
increases with the gas pressure or surface density (Kruijssen 2012). 


Likewise, the emergence of the power-law ICMF emerges natu¬ 
rally from the fractal structure of the ISM (Elmegreen & Falgar- 
one 1996), with a high-mass truncation that is set by the Toomre 
mass of the host galaxy and increases with the pressure (Kim & 
Ostriker 2001; Kruijssen 2014; Adamo et al. 2015). It was shown 
by Krumholz & McKee (2005) that the Toomre mass can be ex¬ 
pressed in terms of simple observables, i.e. the gas surface density, 
angular velocity, and Toomre Q disc stability parameter, if the as¬ 
sumption is made that star formation in galaxies occurs in discs 
residing in hydrostatic equilibrium. 

(iii) While high gas pressures and densities may promote the 
formation of bound stellar clusters by increasing the CEE and the 
ICMF truncation mass, they also effectively destroy bound struc¬ 
ture by tidal perturbations (often referred to as ‘impulsive’ or ‘tidal’ 
shocks) due to encounters with GMCs (e.g. Spitzer 1958, 1987; 
Ostriker, Spitzer & Chevalier 1972; Gieles et al. 2006b). This con¬ 
cept of a high disruption rate at young cluster ages (named the 
‘cruel cradle effect’ by Kruijssen et al. 2011; also see Elmegreen 
& Hunter 2010) is of key importance during the early phase of GC 
evolution within the host galaxy disc, particularly in high-pressure 
environments (e.g. Elmegreen 2010; Kruijssen et al. 2014). Under 
such conditions, other disruption mechanisms such as evaporation 
or dynamical friction destroy stellar clusters more slowly than tidal 
shocks by at least an order of magnitude (Kruijssen et al. 2011, 
2012b). 

(iv) The rapid disruption of stellar clusters by tidal shocks im¬ 
plies that YMCs are unlikely to survive for cosmological time- 
scales unless they migrate out of the disc into the host galaxy halo. 
Present-day GCs are not associated with their host galaxies’ gas 
reservoirs, illustrating that GCs must have gone through a migra¬ 
tion event. ^ This could be caused by several mechanisms, but nu¬ 
merical simulations of the formation and/or disruption of clusters 
during hierarchical galaxy formation suggest that the most obvious 
redistribution of GCs takes place when galaxies merge or are tidally 
stripped by a more massive galaxy (Kravtsov & Gnedin 2005; Pri¬ 
eto & Gnedin 2008; Kruijssen et al. 2012b; Rieder et al. 2013).^ 
These processes ‘liberate’ GCs from the host galaxy disc into the 
halo. The migration time-scale is set by the rate of galaxy mergers 
with similar-mass galaxies (major mergers) or more massive galax¬ 
ies (accretion events), which is /merge = 10^-10^° yr at z ~ 3 
(e.g. Genel et al. 2009) and increases with cosmic time and galaxy 
mass (e.g. Fakhouri, Ma & Boylan-Kolchin 2010). As a result, the 
migration rate of GCs into low-density environments per GC must 
have peaked at high redshift and at low galaxy masses. 

(v) Current numerical simulations of cluster formation and dis¬ 
ruption in a galactic context all find that the survival probabilities 
of pre-existing clusters are indeed greatly increased by the redis¬ 
tribution of matter during galaxy mergers and hierarchical galaxy 
growth (Kruijssen et al. 2012b; Renaud & Gieles 2013; Rieder et al. 
2013). Only the clusters formed at the peak of a merger-induced 
starburst experience greatly enhanced levels of disruption due to 
the high gas densities (Kruijssen et al. 2012b). As such, the surviv¬ 
ing clusters in the merger remnant predate the main starburst and 
typically have ages up to 200 Myr older than the field stars formed 
in a merger (Kruijssen et al. 2011). Note that mergers also play a 


^ Our use of the term ‘migration’ strictly refers to the escape of GCs from 
the part of (6D) position-velocity space that is occupied by the gas. 

^ A more speculative scenario could be that GC migration is facilitated by 
tidal heating or secular evolutionary processes such as dynamical interac¬ 
tions with the gas clumps seen in high-redshift galaxies. 
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key role in setting the properties of the composite GC populations 
observed at z = 0, of which the widely different metallicities and 
binding energies suggest that they formed in a variety of galax¬ 
ies that merged to form the present-day system (e.g. Searle & Zinn 
1978; Forbes & Forte 2001; Muratov & Gnedin 2010; Tonini 2013). 

(vi) Once GCs reside in the host galaxy halo, they slowly dis¬ 
solve by evaporation (Spitzer 1987; Fall & Zhang 2001) at a rate 
that depends on the strength of the tidal field and hence on the 
galactic environment (see Baumgardt & Makino 2003; Gieles & 
Baumgardt 2008, as well as a more extensive discussion of this 
point in §3.5 of Kruijssen 2014). For a flat rotation curve, the dis¬ 
ruption time-scale due to tidal evaporation scales linearly with the 
galactocentric radius (Baumgardt & Makino 2003). As highlighted 
in §2.1.3, evaporation therefore cannot explain the observed insen¬ 
sitivity of the GCMF peak mass and the specific frequency to the 
host galaxy mass or the galactocentric radius (at fixed metallicity), 
unless the strong environmental dependence of evaporation is ig¬ 
nored (as in e.g. Fall & Zhang 2001; Prieto & Gnedin 2008; Katz 
& Ricotti 2014; Li & Gnedin 2014). Secondary disruption mecha¬ 
nisms such as bulge or disc shocks (Gnedin, Hernquist & Ostriker 
1999) are environmentally dependent too. If most GCs are the prod¬ 
ucts of regular YMC formation in high-redshift galaxies, both the 
shape of the GCMF and the specific frequency must therefore have 
been set at early times. As will be shown in §4, the best candidate is 
the rapid disruption phase within the host galaxy disc, i.e. before the 
GCs attained their current spatial distribution within galaxy haloes. 
This conclusion is supported further by the fact that the present- 
day mass loss rates of GCs are too low to facilitate the evolution of 
a Schechter-type ICMF into the peaked present-day GCMF within 
~ 10^° yr (Fall & Zhang 2001; Vesperini et al. 2003; Kruijssen & 
Portegies Zwart 2009). While GC evolution in galaxy haloes can¬ 
not have strongly affected the present-day masses and numbers of 
GCs, it does play a key role in setting the radii of GCs (Gieles, 
Heggie & Zhao 2011; Madrid, Hurley & Sippel 2012). Recent ob¬ 
servational work shows that this structural evolution of GCs within 
galaxy haloes is still ongoing (Baumgardt et al. 2010). 


3 SHAKEN, THEN STIRRED: A TWO-PHASE MODEL 
EOR GLOBULAR CLUSTER EORMATION 

Thus far, models aiming to explain the properties of the GC popula¬ 
tions seen in the local Universe have focused mainly on the dynam¬ 
ical evolution of GCs in the haloes where they presently reside (Fall 
& Rees 1977; Chernoff & Weinberg 1990; Gnedin & Ostriker 1997; 
Vesperini 2001; Fall & Zhang 2001; McLaughlin & Fall 2008; 
Lamers, Baumgardt & Gieles 2010; Muratov & Gnedin 2010; Li 
& Gnedin 2014), thereby omitting the early phase of GC formation 
and evolution within the gas-rich discs of their host galaxies. In this 
section, we use the available observational and theoretical evidence 
discussed in §2.1 and §2.2 to present a simple two-phase model for 
GC formation that includes this early phase of cluster formation 
and evolution as well as the migration of GCs into the galaxy halo 
and their subsequent, quiescent evolution until the present day. 

The model is illustrated graphically in Figure 2 and is ex¬ 
plained fully below. A summary is provided in §3.5, where we also 
highlight the differences with respect to previous models of GC 
formation and evolution in the context of galaxy formation. In §4, 
the model results are compared to observations of present-day GC 
populations, after which the implications of the model as well as its 
testable predictions are discussed in §5. 

The GC formation model presented here connects the forma¬ 


tion of YMCs in the local Universe (§2.1.1), the properties of high- 
redshift, star forming galaxies (§2.1.2), the hierarchical growth of 
galaxies (§2.2), and the physical mechanisms governing cluster 
evolution until the present day (§2.2). The model is constituted by 
the following steps. 


3.1 Initialisation of the host galaxy and its cluster population 

At the time of GC formation, the properties of the host galaxy 
and its cluster population are initialised according to the obser¬ 
vational and theoretical constraints on local YMC formation and 
high-redshift star formation. 


(i) The progenitor clusters of GCs formed at a redshift z, which 
may vary with the galaxy mass and metallicity. For illustration, 
it is here chosen to be z = 3 independently of the galaxy mass, 
corresponding to a look-back time of t = 11.7 Gyr (consistent 
with the ages of GCs, e.g. Wright 2006; Forbes & Bridges 2010; 
VandenBerg et al. 2013). In the context of Figure 2, we thus set 
tform = 0 Gyr and tnow = 11.7 Gyr. This is obviously a simpli¬ 
fication given the extended formation histories of GC populations, 
which is made to keep the model as transparent as possible. In fu¬ 
ture work, we will add the complexity of extended star and cluster 
formation histories (see §6). 

(ii) In the model, the properties of the GC population are de¬ 
rived as a function of the host galaxy dark matter halo mass at GC 
formation Mh. Given some value of Mh, the host galaxy stellar 
mass Mstar follows from the commonly-used abundance matching 
model of Moster, Naab & White (2013). Their abundance match¬ 
ing technique makes the assumption that Mgtar is a monotonically 
increasing function of Mh in order to connect the observed galaxy 
mass function to the halo mass function from the Millennium Sim¬ 
ulation (Springel et al. 2005) for the redshift range z = 0-4. Ana¬ 
lytically, the relation is given by 


Mstar 



( 2 ) 


At the adopted z = 3, the constants are A = 0.033, B = —0.76, 
C = 0.85, and Ml = 3 X 10^^ Mq (Moster, Naab & White 2013). 
Given some stellar mass, the metallicity [Fe/H] of the host galaxy 
(and hence of its stellar cluster population) is obtained through 
equation (1) in §2.1.2. This direct relation to metallicity is impor¬ 
tant, because it shows that the metallicities of GCs trace the masses 
of the galaxies that they formed in. 

(iii) Stellar clusters of metallicity [Fe/H] are generated accord¬ 
ing to a Schechter-type ICMF: 

diV 

— (xM“ exp (-Mi/Me), (3) 

with index a = —2, a. lower mass limit of Mmin = 10^ M© (e.g. 
Lada & Lada 2003; Kruijssen 2014), and a truncation mass Me, 
normalised to a total mass in clusters of FMstar. This is equivalent 
to assuming a roughly instantaneous formation of the host galaxy’s 
stellar and cluster populations. 

(iv) The metallicity of the cluster population and its host galaxy 
is translated to an ICMF truncation mass Me by using the relation 
between the mean metallicity of the present-day GC population and 
the GCMF truncation mass Me (i.e. combining the [Fe/H]-Ms re¬ 
lation from Figure 13 of Peng et al. 2006 with the Me-Ms relation 
from Figure 16 of Jordan et al. 2007). We thus obtain an approxi¬ 
mate relation of 


(c) 2015 RAS, MNRAS 000, 1-29 




8 /. M. D. Kruijssen 


YMC formation 



(a) 


Cluster migration 



(b) 


GC populations 



Time 
-> 


orm Emerge ^now 


1 1 

1 1 

Disc: gas-rich 

Halo: gas-poor 

^disc 

^halo 

Rapid disruption 
by tidal shocks 

Slow disruption 
by evaporation 


Figure 2. This figure shows the time-line of the two-phase "shaken, then stirred' model for GC formation that is presented in §3. At a time tform. stars and 
stellar clusters are formed through the regular star formation process in the gas-rich and high-pressure discs seen at high redshift (see panel a). During their 
subsequent evolution within the host galaxy disc (Phase 1), the clusters undergo rapid disruption by tidal shocks for a total duration fdisc^ until at a time 
f merge the host galaxy undergoes a merger and the clusters migrate into the galaxy halo (see panel b), thereby increasing their long-term survival chances and 
thus becoming GCs. The GCs subsequently undergo quiescent dynamical evolution in the galaxy halo (Phase 2), which is characterised by slow disruption 
by evaporation (see panel c). This phase lasts for a total duration thaio until the present day fnow- Panel (a) shows the clumpy 2 ; = 4.05 galaxy GN20 from 
Hodge et al. (2012), panel (b) shows the migration of stellar clusters in a galaxy merger model from Kruijssen et al. (2012b), and panel (c) shows the GC 
population of the early-type galaxy NGC 4365 from the SLUGGS survey (Brodie et al. 2014). Panels (a) and (c) reproduced by permission of the AAS. 


log (Mc/Mo) - 6.5 + 0.7[Fe/H]. (4) 

Both the [Fe/H]-Ms and Mc-Mb relations are tightly defined, 
hence equation (4) should also have little scatter (<0.3 dex). The 
dependence on galaxy mass (and hence [Fe/H]) of the truncation 
mass must be physical in nature rather than a trivial statistical result 
of lower-mass galaxies producing fewer clusters. When excluding 
nuclear clusters, no dwarf galaxy (Mstar < 5 x 10^ M©) is known 
to host a GC more massive than 10® M© (Georgiev et al. 2010; 
Kruijssen & Cooper 2012), whereas the much more massive Milky 
Way has several such massive GCs without any of the Fe-spreads 
that would indicate a nuclear origin. This shows that even when 
adding up the GC populations of a large number of dwarf galaxies, 
no GCs as massive as in the Milky Way are found, which thereby 
substantiates the Mc-Mb relation of Jordan et al. (2007). 

(v) To obtain the gas surface density and the CFE F at the time 
of GC formation, we use the relation between the truncation mass 
Me and the Toomre mass Mt, i.e. by writing (also see §2.1.1): 

Me = Mei,T = eVMT = er(S, Q, n) , (5) 

where E is the gas surface density, Q is the Toomre (1964) sta¬ 
bility parameter, and Q is the angular velocity. The final equality 


in this expression assumes that young GCs formed in galaxy discs 
in hydrostatic equilibrium to express Mt in terms of E, Q, and Q 
(Krumholz & McKee 2005, also compare to the expression given in 
§2.1.1). We assume a Toomre stability parameter of Q = 1, and use 
the SINS galaxy sample of z = 1.1-3.5 galaxies (Forster Schreiber 
et al. 2009) to obtain a power-law fit of Q as a function of galaxy 
mass Mstar, yielding 

/ A/f \ 0 

“ ="' 

for Mstar = 10^-10^^ ® M©. At fixed Q = 1, this relation implies 
that the mid-plane gas density increases with the host galaxy stellar 
mass as pisM oc (Q/Q)^ oc Ms®tar^ (cf. Krumholz & McKee 2005). 
Finally, we adopt a SFE of e = 0.05 (e.g. Lada & Lada 2003; Evans 
et al. 2009; Kruijssen 2014), and use the CEE model of Kruijssen 
(2012) to obtain F(E, Q, Q). Eor a given galaxy mass Mstar, we 
then numerically solve the implicit relation of equation (5) to obtain 
E (and hence F and Mt). 

(vi) A relation should be adopted between the initial cluster 
masses and their half-mass radii. For simplicity, we adopt a con¬ 
stant radius across the entire cluster population, which is consis¬ 
tent with nearby YMCs (see §2.1.1), and assume that during the 
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initial phase of cluster evolution in the host galaxy disc, the half¬ 
mass radius is proportional to the mean tidal radius of the clus¬ 
ter population, i.e. n oc As shown by equation (6), the 

galaxies under consideration here have a typical angular velocity 
of = 0.1 Myr“^ as opposed to Q = 0.025 Myr“^ in the so¬ 
lar neighbourhood and other nearby spiral galaxies (e.g. Kennicutt 
1998). The characteristic radius of nearby YMCs of rh = 3.75 pc 
(Larsen 2004) thus becomes rn = 1.5 pc in cluster-forming galax¬ 
ies at z ^ 3. 

In summary, we choose a halo mass Mh and GC formation redshift 
(taken to be z = 3) to initialise a galaxy with a known stellar mass 
Mstar, metallicity [Fe/H], angular velocity Q, gas surface density 
E, and mid-plane gas volume density pisM- This galaxy hosts a 
young stellar cluster population with a known ICMF and CFE, and 
we make a physically-motivated assumption for the typical cluster 
radii rn. 


3.2 Shaken by the cruel cradle effect: cluster evolution 
during the rapid-disruption phase 


After the formation of the cluster population, the clusters evolve 
within the host galaxy disc in which they formed (Phase 1). During 
this phase, clusters are disrupted by the tidal injection of energy 
during encounters with GMCs (i.e. tidal shocks) in the gas-rich, 
high-pressure environment. 

We describe the disruption process using the model of Kruijs- 
sen (2012, equation 39 - also see Gieles et al. 2006b), which for 
the assumptions and definitions made in the previous steps (e.g. a 
gas disc in hydrostatic equilibrium) can be expressed in terms of a 
mass loss rate 


/dM\ M 

V tcce’ 

on a disruption time-scale 


(7) 


t 


cce 



PiSM 

M© pc“^ 


f M \ 

VIO® Mo) 




( 8 ) 


where the subscript ‘cce’ refers to the ‘cruel cradle effect’ (see 
§2.2), t 5 ,cce = 176 Myr is the proportionality constant that very 
weakly depends on a cluster’s structural parameters but is taken to 
be constant here, /s = Eqmc/E is the ratio between the GMC 
surface density and the mean gas surface density, and 0ad is the 
‘adiabatic correction’ (Weinberg 1994), which accounts for the ab¬ 
sorption of the tidally injected energy by adiabatic expansion in 
sufficiently dense clusters (see below). Each quantity has been nor¬ 
malised to the typical numbers under consideration here.^ 

Equation (8) shows that clusters are disrupted most rapidly in 
galaxies with high absolute ISM densities as well as high surface 
density contrasts between the GMCs and the diffuse ISM. Because 
the typical disruption time t 5 ,cce is much shorter than a Hubble 
time, the long-term survival of clusters is only possible if they mi¬ 
grate out of the host galaxy disc on a short time-scale. In our model, 
cluster migration is facilitated by galaxy mergers (see below) and 
hence the formation of long-lived GCs requires a high merger rate, 
as is seen at high redshift. Specifically, the long-term survival of 
clusters with masses up to several 10^ M© typically requires that 
migration should take place within ~ 1 Gyr, comparable to the 
typical merger time-scale at z ^ 3. 


^ Note that the linear dependence on cluster mass M is actually one on the 
half-mass density ph. but a constant radius was assumed in §3.1 above. 


Determining /s in equation (8) requires some knowledge 
about the structure of the ISM. By equating the turbulent pressure 
to the hydrostatic pressure in a self-gravitating GMC, Krumholz 
& McKee (2005) derive an expression for the GMC density in 
an equilibrium disc, which in combination with the Toomre mass 
yields the GMC radius and hence its surface density Eqmc • Adopt¬ 
ing a typical GMC virial ratio = 1.3 we obtain 

/e = 3.92 ^ ^ (9) 


where the molecular gas fraction /moi is a function of E as defined 
in equation (73) of Krumholz & McKee (2005), yielding almost an 
entirely molecular ISM (/moi ~ 1) at the surface densities consid¬ 
ered here (see below) and hence the GMC-to-diffuse gas surface 
density ratio is /e ~ 4. 

Einally, the parameter 0ad in equation (8) is the adiabatic cor¬ 
rection, which accounts for the inefficient conversion of the tidally 
injected energy to dynamical mass loss for tidal shocks taking place 
on a time-scale tpert longer than the cluster’s crossing time ter - in 
such a case, (part of) the energy is absorbed by the adiabatic expan¬ 
sion of the cluster. We use the definition of Weinberg (1994) for a 
time-scale ratio <j)t = tpert/tcr and assume the definition of in 
terms of the cluster density ph and the ambient gas density pisM 
from equation (37) of Kruijssen (2012): 


<?!>ad = (1 + 4 h ) 


, „ / Ph/piSM 

V 10^ 


-3/2 


( 10 ) 


Eor the adopted half-mass radius of rn = 1.5 pc and at an 
ISM density of pisM ~ 3 M© pc“^, the adiabatic dampen¬ 
ing of tidally injected energy becomes important (i.e. <j)t > 1) 
for cluster masses M > 10^ M©, which substantially increases 
their survival fractions and may prevent disruption altogether for 
M > several 10^ M©. 

Eor clusters in the initial mass range Mi = 10^-10^ M©, we 
numerically integrate equation (7) to obtain cluster masses M dur¬ 
ing the rapid-disruption phase at some time t after they formed. The 
cluster mass function at time t is obtained by converting the ICME 
of equation (3) through conservation of the number of clusters: 

diV _ diV dMi 

dM “ dMi dM ■ ^ ^ 

In summary, given a cluster mass spectrum and the host galaxy 
ISM properties from §3.1, the cluster masses are evolved by ac¬ 
counting for tidal shocks. Other than a small set of parameters dis¬ 
cussed in Kruijssen (2012) that only weakly affect the mass loss 
rate and depend on e.g. the structural properties of the clusters, this 
step of the model has no free parameters. The main uncertainty 
is the assumption that cluster radii are independent of the cluster 
mass, which is discussed in more detail in §6. 


3.3 Cluster migration into the galaxy halo by galaxy mergers 

The rapid disruption phase continues until the clusters migrate into 
the galaxy halo. The migration agent may be any process that facil¬ 
itates cluster migration, such as secular evolution within the thick, 
gas-rich host galaxy disc, tidal heating, rapidly changing gravita¬ 
tional potentials due to galaxy-wide feedback (Pontzen & Gover- 
nato 2012; Maxwell et al. 2014), or (major or minor) galaxy merg¬ 
ers. In this model, we assume that the migration of clusters into the 
halo results from galaxy mergers. It is shown in §4 that at z ~ 3, 
mergers indeed occur frequently enough to enable the survival of a 
population of young GCs. 
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The merger time-scale of the host galaxy is obtained using the 
galaxy halo merger rates from Genel et al. (2009, see their Figure 3 
in particular), who have used the Millennium Simulation to derive 
the merger rate per unit progenitor halo, time, and galaxy mass ra¬ 
tio as a function of halo mass. We integrate the merger rate over all 
mass ratios x = Mhost/Mother < 3 (where x < 1 indicates that 
the host galaxy is merging with a more massive galaxy). This mass 
ratio range includes major mergers (defined as 1/3 < x < 3), 
as well as the cannibalism of the host galaxy by a more massive 
galaxy (defined as a; < 1/3). Both types of event redistribute clus¬ 
ters into the halo - in the first case (1/3 < x < 3) this is done 
by violent relaxation, whereas in the second case {x < 1/3) the 
redistribution takes place by the tidal stripping of the clusters from 
the host galaxy before it coalesces with the more massive galaxy. 
Mergers with x > 3 are excluded, because these represent cases 
where the GC host galaxy merges with a much less massive galaxy, 
which are unlikely to cause a morphological transformation strong 
enough for the clusters to migrate into the halo (e.g. Hopkins et al. 
2009). 

Inverting the merger rate resulting from the mass ratio integra¬ 
tion yields the merger time-scale, which reflects the cluster migra¬ 
tion time-scale in the model and hence the represents the duration 
of the first, rapid disruption phase tdisc = Emerge — tform (cf. Fig¬ 
ure 2). Genel et al. (2009) provide merger rates at z = 3 for galaxy 
halo masses Mh = { 10 ^^'^, 10 ^^'^, 10 ^^'^}, which are the halo 
masses for which the model results will be shown below (see Ta¬ 
ble 1 in §4.1). 

The only free parameter in this step of the model is the max¬ 
imum halo mass ratio required for cluster migration into the halo 
(taken to be x < 3). As stated above, this choice is well-motivated, 
but we have also verified that the exact choice does not strongly 
affect the merger time-scales. This can be inferred already by vi¬ 
sual inspection of Figure 3(b) of Genel et al. (2009), in which mass 
ratios with a; > 3 represent the clear minority of mergers for halo 
masses Mh < 10^^ M©. The time-scales themselves follow di¬ 
rectly from the Millennium Simulation, and therefore do not repre¬ 
sent free parameters in the context of this model. 


3.4 Stirred by evaporation: cluster evolution during the 
slow-disruption phase 

After their migration into the galaxy halo, the clusters predomi¬ 
nantly lose mass due to evaporation and thus undergo quiescent 
evolution until the present day (Phase 2).^ Following Baumgardt 
(2001) and Lamers et al. (2005), their disruption in the galaxy halo 
can be written as a mass loss rate 


^ The mass loss rate due to evaporation is at least five times higher than 
disc shocking (e.g. Gnedin & Ostriker 1997; Dinescu, Girard & van Altena 
1999; Kruijssen & Mieske 2009; Kruijssen & Portegies Zwart 2009). Bulge 
shocking can be more important and could be modelled by introducing a 
factor 1 — e in equation (13), where e is the orbital eccentricity (Baum¬ 
gardt & Makino 2003). However, it is omitted in our model, because (1) 
there are poor constraints on the orbital eccentricity of GCs as a function of 
metallicity across the galaxy mass range of interest, and (2) for the median 
GC eccentricity in the Galactic halo (e ~ 0.5, see Dinescu, Girard & van 
Altena 1999), the evaporative mass loss rates are only affected by a factor 
of ~ 2. Given that we will show in §4 that the evaporation phase is sub¬ 
dominant for setting the properties of the 2 ; = 0 GC population, omitting 
this aspect is preferable over introducing a free parameter that hardly affects 
the results. 


f 

f^evap 

on a disruption time-scale 

f M 

ievap - is.evap 

where we adopt 7 = 0.7 (e.g. Fukushige & Heggie 2000; Baum¬ 
gardt 2001; Lamers, Gieles & Portegies Zwart 2005; Gieles & 
Baumgardt 2008) and t 5 ,evap oc is the characteristic disrup¬ 
tion time-scale of a M = 10^ M© GC. The disruption time-scales 
of GCs in the Milky Way cover a broad range of t 5 ,evap = 1- 
100 Gyr, depending on the orbital parameters (e.g. Baumgardt & 
Makino 2003; Kruijssen & Mieske 2009). This time-scale is 1-3 
orders of magnitude longer than during the rapid-disruption phase 
in the host galaxy disc (cf. equation 8 ). 

The characteristic time-scale for mass loss by evaporation 
t 5 ,evap is sclectcd such that the near-universal characteristic mass- 
scale of GCs is reproduced as a function of [Fe/H] at z = 0. As 
explained in §4.2 below, this requires t 5 ,evap to decrease with the 
metallicity. For a fiat rotation curve t 5 ,evap oc oc R and thus 
the equivalent assumption is that the GC metallicity [Fe/H] de¬ 
creases with R. As discussed in §2.1.3, GC populations indeed ex¬ 
hibit such a radial metallicity gradient, as metal-poor GCs reside at 
larger galactocentric radii R than metal-rich ones. 

The adopted disruption time-scales can be tested quantita¬ 
tively by assuming a fiat rotation curve and writing 

^5,evap — ^5,0 (77/8.5 kpc), (14) 

with ts,© = 33.8 Gyr for 7 = 0.7 (Kruijssen & Mieske 2009). 
Using the 2010 edition of the Harris (1996) catalogue of Galactic 
GCs, a least-squares fit to the observed gradient of the GC metal¬ 
licity with disruption time (and hence radius through equation 14) 
is 

Fe/H]^bs = -1.11 -0.46 log (t5,evap/10 Gyr). (15) 

In the model, the adopted values of t 5 ,evap for [Fe/H] = 
{—1.1,—0.7,—0.6} are t 5 ,evap/Myr = {15.8,2.5,1.3} (see 
§4.1), implying a least-squares gradient of: 

= -1.03 - 0.50 log (t5,evap/10 Gyr). (16) 

The very similar slopes and normalisations of the observed and 
adopted t5,evap-[Fe/H] relations in equations (15) and (16) show 
that our assumed disruption time-scales give an accurate represen¬ 
tation of the Galactic GC population. Even though t 5 ,evap was cho¬ 
sen to reproduce the observed characteristic GC mass at 2 : = 0, we 
therefore use a relation that is quantitatively consistent with the 
theoretically-predicted disruption time-scales given the observed 
metallicity gradient of Galactic GCs. 

It is plausible that the t5,evap-[Fe/H] (or i 7 -[Fe/H]) gradient 
can be interpreted as an imprint of the binding energy at formation, 
and hence of the galaxy-mass (and hence metallicity) dependent 
redistribution of GCs in mergers and/or tidal stripping at /merge- 
In that case, this relation between [Fe/H] and /5,evap may hold 
universally in galaxies where the decay of GC orbits by dynamical 
friction is negligible (also see §5.9). 

Perhaps surprisingly, an interesting corollary of the above is 
that the large-scale tidal field currently experienced by GCs re¬ 
flects that at earlier times, despite mixing during the further hi¬ 
erarchical assembly of galaxies. After their initial migration, GCs 
may be stripped from their host galaxy if the host galaxy becomes 
a satellite of a more massive halo. However, this does not affect 



/dM\ 

V dij 
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their disruption time-scales, as the stripping takes place when the 
density enclosed by the satellite’s orbit in the more massive halo 
(ptid,sat oc equals the density enclosed by the GC’s or¬ 

bit within the satellite (ptid,GC oc ^^gc)- As a result, the angu¬ 
lar velocity of the GC does not change when it is stripped from 
its host galaxy, nor does stripping alter its disruption time-scale 
f5,evap oc Thc Only exceptions to this rule occur if the satel¬ 
lite is on a very eccentric orbit, which is then inherited by the GC, 
or if the galaxy merger mass ratio is a; ~ 1 (which Genel et al. 2009 
show is rare). The present-day disruption time-scales of GCs may 
thus be taken as a proxy for their disruption time-scales shortly af¬ 
ter their initial migration into the halo (as was shown numerically 
by Rieder et al. 2013). 

After having completed the numerical integration of equa¬ 
tion (7) until the time of cluster migration tmerge, we continue 
evolving the cluster population by integrating equation (12) to ob¬ 
tain cluster masses M during the slow-disruption phase. The clus¬ 
ter mass function at that any time t > tmerge is again obtained by 
converting the ICMF of equation (3) through conservation of the 
number of clusters as in equation (11). 

As indicated in §3. 1 , this evolved cluster population in the host 
galaxy disc still has a single metallicity. However, it is straightfor¬ 
ward to account for the hierarchical growth of galaxies and their GC 
systems by constructing a metallicity-composite population. This 
is done by carrying out a summation of these single-metallicity 
(i.e. single-progenitor galaxy) models with weights given by any 
desired (e.g. observed) distribution function of cluster metallicities. 

In summary, given an evolved cluster mass spectrum at the 
end of the rapid-disruption phase from §3.2, the cluster masses are 
evolved by accounting for tidal evaporation. This step of the model 
is carried out in direct accordance with the results of detailed N- 
body simulations of dissolving clusters. The only free parameter 
(or rather choice) is the relation between metallicity and the char¬ 
acteristic evaporation time-scale. As explained above and in §4.2, 
this relation has been chosen to reproduce the peak mass of the GC 
mass function at z = 0, which we show above is consistent with 
the observed metallicity gradient of GCs in the Milky Way (and 
plausibly in other galaxies, see §5.9). 


3.5 Summary of the model and comparison to previous work 

The simple model presented in this section provides an end-to-end 
description of GC formation and evolution in the context of galaxy 
formation. For a given halo mass, a young stellar cluster population 
forms within a normal, star-forming disc galaxy with known stellar 
mass, metallicity, ISM properties, and cluster population properties 
set by the halo mass, at a fiducial redshift z = 3 (corresponding to 
the median formation redshift of Galactic GCs). The clusters then 
undergo an evolutionary phase characterised by rapid disruption by 
tidal perturbations from molecular clouds and clumps in the host 
galaxy disc. This phase continues until the host galaxy halo merges 
with another halo that has at least 1/3 of the host galaxy’s halo 
mass (i.e. the mass ratio is x = Mhost/Afother < 3), which in 
the model leads to the migration of the clusters into the halo of the 
merger remnant (here assumed to occur instantaneously, either by 
a major merger if x ~ 1 or by tidal stripping if x < 1/3). The 
(progenitor halo) merger time-scale is taken from the Millennium 
Simulation. The further evolution of the cluster population until the 
present day is characterised by slow disruption due to tidal evapo¬ 
ration, with a mass loss rate that is consistent with detailed A^-body 
simulations. 


The total duration of the slow-disruption phase is® thaio = 
tnow — fform “ fdisc, which for GC formation at z ~ 3 (i.e. tnow — 
fform = 11.7 Gyr) and typical merger time-scales fmerge fform — 
foisc < 4 Gyr results in durations at least twice as long as for the 
rapid-disruption phase (thaio/foisc > 2). Despite the fact that GCs 
thus spend most of their lives in the host galaxy halo, the disruption 
rate integrated over this phase is small (compare the characteris¬ 
tic values of t5,cce and t5,evap in §3.2, §3.4, and §4.1 below, with 
f5,cce/f5,evap ^ 0.05). As a rcsult, the long-term survival of clus¬ 
ters becomes much more likely by migrating into the halo. Recall 
that in §1, we defined GCs as a bound stellar cluster that does not 
coincide in phase space with the presently star-forming component 
of the host galaxy. We now see that the migration of GC-like YMCs 
away from the gas-rich host galaxy disc indeed marks the moment 
at which they become long-lived GCs. 

There exist several recent models that connect the origin of 
GCs to galaxy formation (e.g. Prieto & Gnedin 2008; Boley et al. 
2009; Tonini 2013; Katz & Ricotti 2014; Li & Gnedin 2014). How¬ 
ever, the model presented here differs from these previous works by 
introducing two crucial new elements: 

(i) It is the first model that does not make the assumption that 
the (globular) cluster formation rate simply scales with the SFR (or 
that the mass in GCs is proportional to the halo mass), but instead 
introduces an environmentally-dependent fraction of star formation 
occurring in bound stellar clusters, as well as an environmentally- 
dependent ICMF truncation mass. Both of these are in agreement 
with theory and observations of cluster populations in the nearby 
Universe, and the model applies them to the conditions of star for¬ 
mation and the ISM seen in star-forming disc galaxies at high red¬ 
shift. 

(ii) The model is also the first to include the rapid-disruption 
phase in the host galaxy disc due to tidal perturbations from the 
ISM, in a way that is consistent with A^-body simulations of tidally- 
shocked clusters. It will be shown in §4 and §5 that this process 
dominates the total mass loss of the cosmic GC population inte¬ 
grated over a Hubble time. As a result, the evaporation phase is 
sub-dominant, even though it has traditionally been the focus of GC 
formation models and A/^-body simulations of dissolving clusters. 
Effectively, these models thereby immediately place GCs in galaxy 
haloes and skip the (dominant) physics of GC formation and early 
evolution. 

Next to the addition of a detailed model for the formation and 
early evolution of GCs in high-redshift discs, our model also treats 
mass loss by evaporation differently from earlier work (e.g. Gnedin 
& Ostriker 1997; Prieto & Gnedin 2008; Katz & Ricotti 2014; 
Li & Gnedin 2014). These previous studies adopted the expres¬ 
sion from Spitzer (1987), which neglects the non-zero escape time 
of unbound stars (Baumgardt 2001) and assumes that evaporation 
takes place at a rate that is independent of the galactic environment. 
By construction, this has allowed the above previous works to re¬ 
produce the near-universal properties of GC populations (e.g. their 
universal characteristic mass scale, see e.g. Fall & Zhang 2001). 
However, this is inconsistent with numerical simulations of dissolv¬ 
ing clusters (e.g. Baumgardt & Makino 2003; Gieles & Baumgardt 
2008; Lamers, Baumgardt & Gieles 2010), which show that the 
mass loss rate by two-body relaxation is proportional to the square 
root of the tidal field strength (Kruijssen et al. 2011) or inversely 
proportional to the galactocentric radius (for a fiat rotation curve). 


See Figure 2 for the definitions of these time-scales. 
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In summary, the model presented here accounts for the key 
physical mechanisms which are known to be important for cluster 
formation and evolution. As such, it improves on previous work 
in terms of the cluster modelling. However, the coupling to galaxy 
formation in this (largely analytical) model is more simplified, as 
it makes use of a simple merger time-scale and a step function to 
describe the migration of clusters from discs to haloes. While this 
does allow us to easily identify which physical component is re¬ 
sponsible for which aspect of the observables, the obvious benefit 
of previous numerical simulations is that they provide predictions 
for the spatial distribution and kinematics of GC populations. We 
will link our new model to such methods in several future papers. 


4 COMPARISON TO OBSERVED GLOBULAR CLUSTER 
SYSTEMS 

The two-phase model for GC formation presented above is consis¬ 
tent with the formation of YMCs in the local Universe, as well as 
with the star-forming properties of high-redshift galaxies and the 
physics of hierarchical galaxy formation. The obvious next ques¬ 
tion is whether the resulting properties of the GC population agree 
with what is observed. We will now present the results for five dif¬ 
ferent galaxy models, focusing on several of the observables dis¬ 
cussed in §2.1.3. In particular, we will discuss the GC mass func¬ 
tion (§4.2), the specific frequency and metallicity bimodality of 
GCs (§4.3), the relation between the GC population mass and the 
dark matter halo mass (§4.4), and the relation between the number 
of GCs and the supermassive black hole mass (§4.5). 

4.1 Considered example cluster population models 

The above model is used to describe the formation and evolution 
of GCs in three example galaxies with z = 3 halo masses of 
Mh = 10^2.5^ for which Genel et al. (2009) pro- 

vide merger rates. In addition, two models are included to address 
the extremely high mass ratio between GCs and field stars in dwarf 
galaxies like the Fornax dSph galaxy at metallicities [Fe/H] < — 2 
(see e.g. Larsen, Strader & Brodie 2012 and §2.1.3). For these mod¬ 
els, the z = 3 halo mass is taken to be Mh = 10^° M©, corre¬ 
sponding to a stellar mass of 4 x 10® M© (obtained by extrapolating 
the galaxy mass-halo mass relation from Moster, Naab & White 
2013). This is roughly consistent with the total stellar mass of For¬ 
nax at metallicities [Fe/H] < — 2 (Larsen, Strader & Brodie 2012). 
The mean metallicity is taken to be [Fe/H] = —2.3, and we adopt 
a ICMF truncation mass of Me = 7 x 10® M©. The disruption 
time-scale during the halo phase is chosen to be t 5 ,evap ~ W yr, 
consistent with a galactocentric radius of R — 2 kpc (cf. Mackey 
& Gilmore 2003) in the observed gravitational potential (Walker & 
Penarrubia 2011). 

The two Fornax models differ only in one parameter. The 
merger time-scale of the first Fornax model is obtained by extrap¬ 
olating a power-law fit to the Genel et al. (2009) relation between 
Emerge and Mh down to a halo mass of Mh = 10^® M©, resulting 
in tmerge = 1.4 Gyr. The merger time-scale of the second Fornax 
model is chosen to reproduce the high observed GC-to-field star 
mass ratio by limiting the duration of the rapid-disruption phase in 
the host galaxy disc, resulting in tmerge = 0.1 Gyr. This corre¬ 
sponds to a couple of orbital times. In the context of our model, 
such a short time-scale represents a scenario in which the high spe¬ 
cific frequency at low metallicities in Fornax may have been caused 
by e.g. the formation of the GCs during the first pericentre passage 


of the merger that eventually triggered their migration into the halo 
(see below). 

Table 1 shows the defined properties (indicated by an asterisk) 
of the galaxy models discussed in this section, as well as those de¬ 
rived using the physical descriptions or fitting functions described 
above. Having initialised the cluster population in a way consis¬ 
tent with observations of stellar cluster formation in the local Uni¬ 
verse, the key question now is whether the derived galaxy prop¬ 
erties are consistent with the ‘normal’, star-forming disc galaxies 
observed at redshift z > 2 (e.g. Tacconi et al. 2010). The gas sur¬ 
face densities in Table 1 place these galaxies at the low end of 
the surface density range of gas-rich, normal galaxies at z > 2 
(Tacconi et al. 2013), indicating that the conditions needed for GC 
formation are indeed common at high redshift. The correspond¬ 
ing pressures are also modest, falling in the range Pjk — 10®- 
10^ K cm“®. Although this is 2-3 orders of magnitude higher than 
the turbulent and hydrostatic pressures in the Milky Way disc, it is 
within the range of pressures observed at > 2, which reach max¬ 
ima of Pjk — 10^-10® K cm“® (see §2.1.2 and e.g. Swinbank 
et al. 2012). The resulting gas volume densities are up to 100 times 
higher than in the Milky Way disc, but again consistent with the 
prevalent high-density gas tracer emission detected in high-redshift 
galaxies (e.g. from high-J CO transitions, which suggest the pres¬ 
ence of large mass reservoirs with peak densities p > 50 M© pc“®, 
see Danielson et al. 2013). As Table 1 shows, the increase of the gas 
pressure is accompanied by a corresponding increase of the fraction 
of star formation occurring in bound clusters F, as well as a strong 
decrease of the cluster disruption time-scale tece- 

Of course, observed z = 0 galaxies will host a composite GC 
population with a range of metallicities and formation redshifts, 
rather than the ‘Simple GC Populations’ of Table 1, which each 
have a single age and metallicity. The galaxy models presented here 
are chosen to cover the GC metallicity range. If GCs are indeed the 
products of regular YMC formation at high redshift, these galaxy 
models should reproduce the observed properties of the present- 
day GC population as a function of metallicity. In the model, we 
account for the further growth and composite nature of GC pop¬ 
ulations by carrying out a summation over a set of models with 
different metallicities (see §3.4 and §4.2). 

4.2 Two-phase cluster disruption and the globular cluster 
mass function 

The left-hand panel of Figure 3 shows the shape of the cluster mass 
function for the five galaxy models of Table 1 at the three key mo¬ 
ments of cluster formation, their migration, and the present day 
(tform, Emerge, and tnow as defined in Figure 2 - note that tmerge 
is different for each of the models because it increases with halo 
mass and metallicity). The figure shows that the Schechter-type 
ICMFs with slightly different truncation masses Me evolve to a 
peaked shape during the rapid-disruption phase in the host galaxy 
disc, due to the tidal disruption of low-mass clusters. The peak 
mass varies substantially with the host galaxy, but is almost univer¬ 
sally larger than the present day peak mass (i.e. Mpeak(Emerge) > 
Mpeak(tnow) = 2 X 10® M©). The possible increase of the peak 
mass to even higher masses is arrested by the decrease of the 
dynamical time with cluster density (and with mass, given that 
Th oc M^ with 13 < 1/3, see §2.1.1), which for clusters beyond 
the peak mass becomes so short that they can respond adiabatically 
to tidal perturbations and do not suffer the violent mass loss that de¬ 
stroys the lower-mass clusters. This transition from impulsive per¬ 
turbations to adiabatic expansion causes a slight steepening of the 
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Table 1. Defined (*) and derived properties of the galaxy models at 2 ; = 3 



Fornax 

[Fe/H] = -1.1 

[Fe/H] = -0.7 

[Fe/H] = -0.6 

log(Mh/M 0 ) 

10 . 0 * 

11.3* 

12.5* 

13.5* 

log (Mstar/M©) 

6.6 

8.9 

10.7 

11.1 

[Fe/H] 

-2.3* 

- 1.1 

-0.7 

- 0.6 

Emerge [Gyr] 

{1.4,0.!*} 

2.0 

2.3 

3.4 

log(Mc/M©) 

5.8* 

5.7 

6.0 

6.1 

log (Mt/M©) 

7.6 

7.3 

7.6 

7.6 

P/k [10® K cm-®] 

1.06 

2.09 

6.54 

8.78 

E [M© pc-2] 

103 

145 

257 

297 

n [Myr-i] 

0.06 

0.09 

0.12 

0.13 

PISM [M© pc-®] 

0.82 

1.75 

3.21 

3.70 

/e 

4.10 

4.01 

3.95 

3.94 

r 

0.38 

0.46 

0.55 

0.56 

tcce(10® M©) [Gyr] 

2.52 

0.36 

0.09 

0.06 

Gvap( 10 ® M©) [Gyr] 

94.87 

15.81* 

2.53* 

1.26* 


dashed curves at M ~ 10^ M©. In addition, the further increase of 
the peak mass due to disruption is also decelerated because it runs 
into the exponential truncation of the ICMF (cf. Gieles 2009). 

During the subsequent slow-disruption phase in the galaxy 
halo, the peak mass decreases due to the gradual, evaporation- 
driven mass loss of the clusters that survived long enough to mi¬ 
grate into the halo. As discussed in §3.4, we have assumed a mass 
loss rate due to evaporation such that the observed Mpeak(tnow) ~ 
2 X 10^ M© is reproduced. At a time tmerge (dashed lines in Fig¬ 
ure 3), higher-metallicity models have higher peak masses, hence 
our assumption implies that high-metallicity GCs must lose more 
mass due to evaporation, i.e. the corresponding disruption time- 
scale t 5 ,evap must dccreasc with metallicity to move the peak 
masses back to a similar mass-scale by z = 0. We discussed in §3.4 
that the resulting adopted relation of a decreasing t 5 ,evap oc oc 
R with [Fe/H] is quantitatively consistent with the observed metal¬ 
licity gradient of the Galactic GC population. This implies that even 
though t 5 ,evap was tuucd to reproduce the observed Mpeak, the re¬ 
sult is equivalent to adopting the theoretically predicted evaporation 
time-scales. This is an important test of the model. 

The above results show that our model predicts a near- 
universal peak mass of the GCMF at z = 0. The rapid destruc¬ 
tion of low-mass clusters in their host galaxy discs is most effective 
at high galaxy masses and metallicities, causing the universally- 
high Mpeak to increase with metallicity by the end of the rapid- 
disruption phase (i.e. at tmerge). This Mpeak-[Fe/H] relation is 
subsequently washed out by the opposite metallicity-dependence 
of mass loss due to evaporation in the galaxy halo, in which metal- 
rich GCs reside at smaller galactocentric radii and lose more mass 
than metal-poor GCs, resulting in a near-universal peak mass af¬ 
ter a Hubble time. As discussed in §5.9, the assumption under¬ 
pinning this prediction is that the relation between [Fe/H] and 
t5,evap oc OC R obscrvcd in the Milky Way also holds in 

other galaxies to within a factor of several. This should be expected 
if the relation reflects a more fundamental relation between the GC 
metallicity and the binding energy within the host galaxy at the time 
of GC formation, which in turn would simply reflect the galaxy 
mass-metallicity relation. 

In our model, the GC population is divided into metallicity 
bins to enable a direct connection to the properties of the host 
galaxy at the time of GC formation through the galaxy mass- 
metallicity relation. However, the GC populations of observed 
galaxies are constituted by a broad range of metallicities (see 


§2.1.3). The right-hand panel of Figure 3 shows the evolution from 
the Schechter-type ICMF at z = 3 to the peaked GCMF at z = 0 
for a composite-metallicity GC population representative for the 
Milky Way. This panel accounts for the hierarchical growth of the 
metallicity-composite Galactic GC population by combining the 
GCMFs for the different metallicities shown in the left-hand panel, 
each scaled by the number of Galactic GCs in the corresponding 
metallicity bin (using the 2010 edition of the Harris 1996 catalogue, 
with bin separations at [Fe/H] = {—1.7, —0.9, —0.65}). The flg- 
ure shows that at z = 3, our model matches the ICMF of YMCs 
observed in the Antennae galaxies, whereas it reproduces both the 
normalisation and shape of the GCMF of the Milky Way at z = 0 
(note that this model has not been fit to reproduce the data).^ The 
near-universal peak mass was already evident when considering the 
different metallicities individually, but for composite-metallicity 
GC populations any remaining, small variations average out and 
the universality of the peak mass becomes almost inescapable. 

The two-phase evolution of the GCMF contrasts strongly with 
the classical idea that the evolution from the power-law ICMF to a 
peaked GCMF is caused by evaporation only. In that model, haloes 
with low GC evaporation rates should have low peak masses. In the 
two-phase model, there exist no GC systems with low peak masses 
because most of the evolution towards a peaked GCMF took place 
in the first few Gyr after the clusters were formed. In our model, 
the approximate universality of the peak mass disappears only at 
metallicities where the present-day peak mass is not attained during 
the rapid-disruption phase [i.e. Mpeak (Emerge) < 2 x 10^ M©]. We 
And this occurs at metallicities of [Fe/H] ^ —1.5. Towards lower 
metallicities (and hence galaxy masses), we should thus expect a 
trend of decreasing Mpeak- This is remarkably consistent with the 
observed, weak trend between Mpeak and GC [Fe/H], which sets 
in at Mb ^ —17 and hence [Fe/H] < —1.3 (as seen by combining 
Figure 13 of Peng et al. 2006 with Figure 15 of Jordan et al. 2007). 

As shown by the thick blue solid line, the two-phase model 
also allows for the existence of Fornax-like galaxies with a high 
number of GCs per unit stellar mass at [Fe/H] < —2, provided 
that tdisc = Emerge “ tform was short. In this model, the host galaxy 

^ The slight discrepancy at M > 10® M© may exist because the model 
does not include the former nuclear clusters of tidally stripped dwarf galax¬ 
ies (see § 6 ), of which there may be a handful of examples in the Galactic 
halo (such as ccCen, see e.g. Lee et al. 2009). However, the difference is 
also consistent with the uncertainties on the model parameters. 
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Figure 3. This figure shows that the combination of regular cluster formation at 2; > 2 with subsequent, two-phase disruption yields a peaked GCMF with a 
near-universal peak mass, consistent with observations of present-day GC populations. Left panel: Shown is the number of GCs per unit log M for the five 
different galaxies from Table 1. The legend indicates the colour coding, where the {thin,thick} blue lines represent the {long,short}-fmerge Fornax model 
(towards the bottom right and top left in the figure, respectively). Dotted lines show the initial cluster mass function at formation (fform)^ dashed lines show 
the cluster mass function at the end of the rapid-disruption phase (tmerge), and solid lines show the GC mass function at 2; = 0 (fnow)- The peak mass reaches 
a maximum at fmerge due to the destruction of low-mass clusters in the host galaxy disc, after which it decreases due to the further evaporation of the massive, 
surviving GCs. Right panel: Shown is the number of GCs per unit log M for a composite-metallicity GC population representative for the Milky Way. The 
black curves are obtained by a linear combination of the single-metallicity and single-age GC populations shown by the coloured lines in the left-hand panel, 
scaling each single-metallicity GCMF by the number of Galactic GCs in the corresponding metallicity bins (see the text). At the lowest metallicities, the mean 
GCMF of the two Fornax models was used. The filled dots with Poisson error bars show the observed GCMF of the Milky Way (Harris 1996, 2010 edition), 
whereas open dots show the ICMF of YMCs in the Antennae galaxies (Zhang & Fall 1999, multiplied by a factor of 30). Both provide a good match to the 
modelled cluster mass functions, illustrating how observed GCMFs can be reproduced by a linear combination of the single-metallicity model GCMFs. 


experienced a merger some 10^ yr after the formation of the young 
GCs. Because the densest and most massive clusters can respond 
adiabatically to tidal perturbations, only clusters below a critical 
density or mass (as rn oc with [3 < 1/3, see §2.1.1) get de¬ 
stroyed, whereas clusters above that critical mass hardly lose any 
stars. For our short-tmerge Fornax model, this transition occurs at 
M ~ 10^-® Mq. 

While a short migration time-scale ensures the survival of a 
larger number of GCs than for normal migration time-scales, the 
shape of the observed GCMF in Fornax is not fully consistent with 
that shown here. The observed peak mass is higher than in the 
model by a factor of ~ 3. This could potentially be caused by 
dynamical friction - if the present-day orbits of the GCs (which 
were used to determine t 5 ,evap) were wider in the past, the total 
amount of mass loss during the halo phase may have been minimal. 
However, the low peak mass is already present at tmerge, suggest¬ 
ing that the discrepancy already arises during the rapid-disruption 
phase. It could be that we underestimated the effect of disruption. 
Alternatively, the short-tmerge model was designed to match the 
high observed specific frequency of Fornax (see below), but did so 
while keeping the minimum cluster mass constant. If the minimum 
cluster mass was higher in Fornax, the observed specific frequency 
would require more cluster disruption (i.e. a later migration of the 
GCs into the halo), which would increase the peak mass. 

Finally, we note that the difference between the GCMFs of the 
two Fornax models illustrates the maximum impact of the stochas- 
ticity of the merger process. The adopted time-scale for the long- 
tmerge modcl (thin blue lines in Figure 3) represents a cosmic av¬ 


erage, but since galaxy merging is a Poisson process, deviations to 
much shorter (thick blue lines in Figure 3) or longer merger time- 
scales (and hence tdisc) are possible. The importance of Poisson 
noise should decrease with galaxy mass, because the GC popula¬ 
tions of more massive galaxies should originate from a larger vari¬ 
ety of progenitor galaxies. 

4.3 The specific frequency as a function of metallicity and 
stellar mass 

Figure 4 shows the number of GCs per unit galaxy stellar mass [the 
specific frequency Tn = iVGc(Mstar/10® M©) as a function 
of the host galaxy metallicity and stellar mass for the five galaxy 
models of Table 1 at the time of GC migration and at the present 
day (tmerge and tnow, see Figure 2). The figure shows that the 
specific frequency decreases with the metallicity and host galaxy 
mass, which is caused by an increase of both the disruption rate and 
merger time-scale (i.e. the duration of the rapid-disruption phase) 
with galaxy mass (see Table 1). We find that this decrease is quan¬ 
titatively consistent with the observed trend (also see below). 

In §2.1.3, we have drawn particular attention to the relation 
between Tn and the host galaxy mass at the time of GC formation 
tform, because its observed invariance with galactocentric radius 
implies that it cannot be affected by GC evaporation in the galaxy 
halo. It must therefore have been set shortly after GC formation. 
Figure 4 shows that this is indeed the case - the specific frequencies 
are in place at tmerge and evolve very little during the halo phase un¬ 
til tnow This is consistent with the result shown in Figure 3 above. 
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Figure 4. This figure shows that the GC specific frequency decreases with 
host galaxy mass and metallicity, because the destruction of clusters dur¬ 
ing the rapid-disruption phase in their host galaxy discs both proceeds more 
rapidly and lasts longer in massive galaxies. Shown is the specific frequency 
Tiv = A^gc (Afstar/lO® Mq ) “ ^ as a function of the host galaxy metallic¬ 
ity (bottom x-axis) and stellar mass (top x-axis), for the five different galax¬ 
ies from Table 1. The legend indicates the colour coding of the coloured 
symbols, where the {open,filled} blue symbols represent the {long,short}- 
f merge Fomax model. Coloured circles show the specific frequencies at the 
end of the rapid-disruption phase (tmerge), whereas coloured squares show 
the specific frequencies at 2 : = 0 (fnow)- The grey line indicates the ana¬ 
lytically predicted relation Tn oc (see the text). The figure shows 

that the relation of decreasing T/v with [Fe/H] and Mstar is established 
at fmerge due to rapid cluster disruption, after which the number of GCs 
remains roughly constant except for the most metal-rich GCs. The lack of 
evolution during the halo phase is consistent with the conclusion of §2.1.3 
that the specific frequency of GCs does not depend on the galactocentric 
radius and hence must have been set during their formation or early evolu¬ 
tion. For comparison, we again include the GC populations of NGC 5128 
and Fornax as black symbols (see Figure 1). 

where we discussed how the number of clusters decreases by tidal 
destruction during the rapid-disruption phase, whereas during the 
halo phase the masses of the surviving GCs decrease and few GCs 
are being destroyed altogether. Consequently, Ngc and Tn are set 
at tmerge (i-c. at z = 1-2), whcrcas the GC masses and radii un¬ 
dergo some mild further evolution within the host galaxy haloes. 

It can be shown analytically how the relation between the spe¬ 
cific frequency and the host galaxy mass and metallicity is shaped 
by cluster disruption. The specific frequency represents the num¬ 
ber of surviving GCs per unit stellar mass and hence combines a 
large part of all the cluster formation and disruption processes dis¬ 
cussed in §2.2. However, the number of YMCs (i.e. young GCs) 
initially formed per unit stellar mass shows little variation with 
galaxy mass (see Table 1), because it only depends on (the weakly 
galaxy mass-dependent) F and Me. The decrease of the specific 
frequency with galaxy mass must therefore refiect an increasing 
trend of the amount of cluster disruption (i.e. a decreasing survival 
fraction). Table 1 shows that the rate and duration of cluster disrup¬ 
tion during the rapid-disruption phase in the host galaxy disc does 
indeed exhibit a strong dependence on the host galaxy mass. The 
disruption time-scale due to tidal shocks tece is expressed in equa¬ 
tion (8), which refiects how rapidly clusters are being disrupted, 
whereas the migration (i.e. merger) time-scale Emerge refiects how 
long clusters are being disrupted. The ratio between these dis¬ 
ruption and merger time-scales thus probes the survival fraction. 


i.e. Tn oc tece/tmerge- The tidal disruption time-scale scales with 
the host galaxy’s stellar mass at z = 3 as tece oc oc Mg{ar^^ 

(see §3.1), whereas the merger time-scales of Table 1 satisfy a 
power-law relation of tmerge = 2.4 Gyr (Mstar/10^° 

As a result, the specific frequency should decrease with stellar mass 
as Tn oc which is indeed seen in Figure 4. 

The above line of reasoning shows that the fundamental rela¬ 
tion of Figure 4 is one between Tn and the host galaxy mass at the 
time of GC formation, because both the gas density and the merger 
rate primarily depend on the galaxy mass in our model. However, 
the information on the initial host galaxy mass is lost during the 
subsequent growth and merging of the galaxies hosting GCs. To 
still trace the mass-dependence of the specific frequency in the lo¬ 
cal Universe, the galaxy mass-metallicity relation at the formation 
redshift of GCs can be used to translate the Tn-M star relation to a 
correlation between Tn and [Fe/H]. This explains the origin of the 
observed TN-[Te/B] relation. 

For metallicity-composite GC populations, the monotonic de¬ 
crease of Tn with galaxy mass and metallicity shown in Figure 4 
should be reversed at high galaxy masses through hierarchical 
galaxy growth. Massive, metal-rich galaxies accrete a large frac¬ 
tion of their mass (20-50% for Mh = lO^^-lO^"^ M©, see Moster, 
Naab & White 2013) in the form of low-metallicity galaxies with 
high Tn, implying that the specific frequency of the massive galaxy 
must increase. Because the accretion rate of low-mass galaxies in¬ 
creases with galaxy mass, this can lead to a minimum of the z = 0 
relation between the specific frequency and the stellar mass. The 
galaxy mass at which this minimum is reached corresponds to the 
point where the GC population becomes dominated by the metal- 
poor GCs accreted from stripped dwarf galaxies. Figure 8 of Peng 
et al. (2008) shows that this is the case for Mb < —22, which in¬ 
deed roughly corresponds to the upturn of the specific frequency 
towards high galaxy masses (Peng et al. 2008, top-left panel of 
Figure 6). For such a composite population of GCs with different 
metallicities (and possibly formation redshifts), the galaxy-wide 
specific frequency can be obtained by integrating the TAr-[Fe/H] 
relation of Figure 4 over the halo field star metallicity distribution 
function. This integration should account for the halo star formation 
history and the dependence of Tn on the formation redshift. For 
instance, the survival fractions of newly-formed clusters decrease 
strongly with cosmic time due to the decreasing galaxy merger rates 
and Toomre masses, whereas the stellar mass may continue to grow 
by in-situ star formation. 

Analogously to the GCMF in §4.2, the Tn difference be¬ 
tween the two Fornax models illustrates the maximum impact of the 
stochasticity of the merger process. While the long-tmerge model 
yields a specific frequency of Tat ~ 200, the short-tmerge model 
reproduces the observed specific frequency of Tn ~ 2 x 10^. This 
order-of-magnitude difference agrees well with the scatter of the 
specific frequency at low galaxy masses (e.g. Peng et al. 2008; 
Georgiev et al. 2010) and should decrease with galaxy mass due 
to the increasing importance of (more numerous and hence less 
stochastic) minor mergers. Despite the model’s good agreement 
with the Fornax dSph, Figure 3 shows that the mean GC mass 
in the model is a factor of 2-3 lower than observed in Fornax 
(M ~ 10^ M©). As a result, the number of GCs per unit stellar 
mass is well-reproduced by the model, but GCs constitute only 9% 
of all coeval stellar mass at z = 0, about a factor of two lower than 
the observed percentage (Larsen, Strader & Brodie 2012). None the 
less, the model does quite well in describing the statistics of For¬ 
nax’ GC population, showing that a short migration time-scale is 
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a feasible hypothesis for explaining the high specific frequency at 
[Fe/H] < — 2 in the Fornax dSph. 

In §2.1.3, it was posited that the Tat -Mstar(tform) relation 
of Figure 4 is of fundamental importance to a theory of GC for¬ 
mation, because it implies several of the other characteristics of 
GC populations. For instance, Tonini (2013) used the z = 0 re¬ 
lation between Tn and Mstar (which is a metallicity and forma¬ 
tion redshift composite of the specific frequencies shown in Fig¬ 
ure 4) to initialise GC populations at redshifts z = 2-4 and show 
that hierarchical galaxy formation naturally produces the observed 
colour (and hence metallicity) bimodality observed in the GC pop¬ 
ulations of massive (Mstar ^ 10® M©) galaxies (see §2.1.3). 
The Tonini (2013) model works, because the specific frequency in¬ 
creases steeply towards low galaxy masses and metallicities, which 
means that accreted dwarf galaxies can boost the number of low- 
[Fe/H] GCs even if their impact on the field star metallicity dis¬ 
tribution is negligible. However, this argument relies on the ad- 
hoc assumption that the Tat -Mstar relation was already in place 
at z ~ 2. Figure 4 shows that this is indeed the case, because the 
GC destruction that sets the Tat -Mstar relation takes place before 
the GCs have migrated into the halo. Combining the physical model 
for the origin of the Tat -Mstar relation presented here with the re¬ 
sult from Tonini (2013), we have thus captured the physics behind 
the bimodality of GC colour and metallicity observed at = 0. 

4.4 The globular cluster system mass per unit halo mass 

Figure 5 shows the relation between the dark matter halo mass Mh 
and the total mass of the galaxy’s GC population MGc,tot for the 
five galaxy models of Table 1 at the time of GC migration and at the 
present day (tmerge and tnow, see Figure 2). As discussed in §2.1.3, 
this relation has received considerable attention in the recent ob¬ 
servational literature. The near-linear slope (i.e. the constant ratio 
log r] = MGc,tot/Mh = —4.5, Durrell et al. 2014) is quite sugges¬ 
tive of a common origin of GCs and dark matter haloes (e.g. Spitler 
& Forbes 2009). However, Figure 5 shows that the model presented 
here also gives rise to a constant GC system mass per unit halo 
mass, even though the GC formation rate is not assumed to scale 
with the host halo mass. It will be shown below that this relation 
follows naturally from the combination of GC disruption (through 
the Tat-M star (tform) rclatiou from Figure 4) and the physics of 
galaxy formation within dark matter haloes. 

The modelled MGc,tot-Mh relation of Figure 5 agrees well 
with the observed relations (Spitler & Forbes 2009; Hudson, Har¬ 
ris & Harris 2014; Durrell et al. 2014), with a normalisation that is 
similar to within a factor of two. The only discrepancy exists for the 
[Fe/H] = —0.6 galaxy, but this should disappear in the composite 
GC populations of z = 0 galaxies, because they are constituted by 
a mix of different metallicities, and therefore contain only a small 
fraction of [Fe/H] > —0.6 GCs. The difference between the latest 
observations (Hudson, Harris & Harris 2014; Durrell et al. 2014) 
and the earlier relation from Spitler & Forbes (2009) is caused by 
a difference in spatial coverage. The Next Generation Virgo Clus¬ 
ter Survey (Ferrarese et al. 2012) as well as the tensing results from 
Hudson, Harris & Harris (2014) now cover a large sample of galax¬ 
ies all the way out to their virial radii, implying a slight increase of 
the dark matter mass per unit GC mass relative to Spitler & Forbes 
(2009). This upward shift is in agreement with the model results. 

We can construct a simple metallicity-composite GC popu¬ 
lation by assuming that ~ 1/6 of the GCs has [Fe/H] > —0.6 
(e.g. Harris 1996; Peng et al. 2006). This results in predicted val¬ 
ues of log ?7 = —4.2 at tmerge and log ?7 = —4.6 at tnow, which 



^GC.tot 

Figure 5. This figure shows that the total GC system mass per unit host 
galaxy halo mass is constant. Shown is the dark matter halo mass as 
a function of the total mass of the GC population, for the five differ¬ 
ent galaxies from Table 1. The legend indicates the colour coding, where 
the {open,filled} blue symbols represent the {long,short}-tmerge Fornax 
model. Coloured circles show the halo and total GC masses at the end of the 
rapid-disruption phase (tmerge), whereas coloured squares show the halo 
and total GC masses at 2 ; = 0 (tnow)- The solid line shows the observed re¬ 
lation from the Next Generation Virgo Cluster Survey (Durrell et al. 2014), 
which covers galaxy haloes out to the virial radius. Other observational es¬ 
timates are shown as the dashed (Hudson, Harris & Harris 2014) and dotted 
(Spitler & Forbes 2009) lines. The figure shows that there is little evolution 
of 77 = MGc,tot/Mh between tmerge and tnow, consistent with the sim¬ 
ilar lack of evolution shown in Figures 3 and 4. Only for the most metal- 
rich GCs there is a deviation, but in observed 2 ; = 0 galaxies this effect 
is washed out because most GCs have metallicities [Fe/H] < —0.6, for 
which this figure shows that 77 constant. As is discussed in the text, the 
constant GC system mass per unit halo mass arises because the decrease of 
the specific frequency with galaxy mass (see Figure 4) is cancelled by the 
increase of the stellar mass fraction with galaxy mass. 

are both within a factor of 2-3 of all observed GC population-to- 
halo mass relations. This is remarkable for two reasons. Firstly, 
we have not assumed any intrinsic connection between GCs and 
dark matter haloes (such as the idea that GC formation took place 
during reionization as the first haloes collapsed). In the model, 
GCs are the products of regular star and cluster formation in high- 
redshift galaxies, yet a one-to-one relation between GCs and haloes 
emerges naturally. Secondly, the model is very simple and omits 
e.g. the composite formation redshifts of 2 : = 0 GC populations and 
further galaxy growth between tmerge and tnow- The good agree¬ 
ment between theory and observations suggests that any mix be¬ 
tween GC sub-populations at different ages and metallicities should 
retain the constant value of 77 . Indeed, the above example of a 
metallicity-composite GC population illustrates how hierarchical 
galaxy growth will further decrease the scatter of 77 across the 
galaxy population due to the central limit theorem. 

The obvious question is what causes the ratio MGc,tot/Mh 
to be constant. The model allows us to identify the reason by ex¬ 
pressing 77 in terms of the TN-Mstar relation seen in Figure 4. By 
writing 

_ -MQCjtot -VgC -VgC -Mstar rp -Mstar /i'7\ 

we see that 77 is the product of Tn oc (see §4.3) and the 

stellar-to-halo mass ratio. For Mh < 10^^'^ M©, Moster, Naab & 
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White (2013) show that the relation between galaxy stellar mass 
and halo mass at z = 3 is Mh oc resulting in Mgtar/Mh oc 

Mstar- Substituting this into equation (17), we indeed find that 77 is 
nearly independent of the galaxy mass with 77 oc 

Present-day galaxies are constituted by composite stellar pop¬ 
ulations assembled through the merging of the single-age and 
single-metallicity populations considered here. Between z = 3 and 
z = 0 , this merging process must have further diluted the already- 
weak correlation between 77 and galaxy mass seen in the single- 
metallicity models of Figure 5. For halo masses Mh > 10^^'^ M©, 
the stellar mass fraction decreases with the halo mass, in stark con¬ 
trast with the increase at lower masses that leads to a constant value 
of 77 through equation (17), indicating that the GC system mass per 
unit halo mass should decrease in massive haloes. However, the GC 
populations in massive haloes are generally observed to have mean 
metallicities [Fe/H] < —0.7 (Peng et al. 2006), indicating that the 
majority of their GCs originates from lower-mass galaxies. As a re¬ 
sult, the ratio between the GC population mass and the halo mass 
is still dominated by GCs that formed in haloes with Mh < 10^^'^ 
for which Mgtar/Mh oc Mg^^r and hence 77 ~ constant. Only for 
galaxy masses of Mstar > 10^^'^ M© at z = 0 (corresponding 
to Mh > 10^^'^ M©), the mean metallicity of the GC population 
reaches [Fe/H] ~ —0.7, possibly leading to a slight steepening of 
the Mh-MGc,tot relation. There is a tentative hint for such a steep¬ 
ening in Figure 3 of Spitler & Forbes (2009), but a larger sample 
of high-mass galaxies is necessary to confirm this. Such a sample 
should also enable the identification of second-order dependences, 
such as on GC sub-population metallicity and galaxy morphology 
(Harris, Harris & Hudson 2015), which may reflect differences in 
the formation histories of the host galaxy and its progenitors. 

In summary, our model suggests that a constant GC system 
mass per unit halo mass arises naturally because both the number 
of surviving GCs per unit stellar mass (Figure 4) and the halo mass 
per unit stellar mass (Moster, Naab & White 2013) decrease with 
the host galaxy stellar mass at the same rate for halo masses Mh < 
10^^'^ M©. When taking the ratio between MGc,tot and Mh, the 
dependence on the stellar mass cancels, resulting in 77 ~ constant. 
This shows that the observed scaling is a ‘coincidence’ and can be 
understood without requiring special, early-Universe physics. 

4.5 The scaling between the number of globular clusters and 
the supermassive black hole mass 

Using a sample of nearby early-type galaxies, it was recently shown 
by Burkert & Tremaine (2010) that the number of GCs scales 
with the mass of the host galaxy’s central supermassive black hole 
(SMBH) M„ i.e. Nac oc Follow-up work by Harris 

& Harris (2011) has confirmed this result, although observations of 
SMBHs in spiral galaxies suggest that the relation may be signif¬ 
icantly shallower in late-type systems (Lomeli et al. in prep.). Ir¬ 
respective of the exact slope of the A^gc-M, relation, it is almost 
certainly not causal in nature - there is no known physical mecha¬ 
nism that directly connects a galaxy’s GC population to the central 
SMBH, other than the suggestion that both classes of objects as¬ 
sembled mainly at early cosmic times (e.g. Kruijssen 2014; Sijacki 
et al. 2015). None the less, we can test if our model reproduces 
the empirical A^gc-M* relation by combining the observed scaling 
relations between the SMBH mass and the host galaxy properties 
with the predicted relations between the number of GCs and the 
host galaxy properties. 

The SMBH mass is primarily related to the velocity disper¬ 
sion or stellar mass of the host galaxy spheroid (Magorrian et al. 


1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000), i.e. a spe¬ 
cific galactic component. Unfortunately, relations specific to a sin¬ 
gle component are not suitable to test whether the A^gc-M, re¬ 
lation is consistent with the GC formation model presented here, 
because our model does not distinguish between the morphological 
components of the host galaxy. Instead, we can use the relation be¬ 
tween the dark matter halo mass Mh and the SMBH mass of Mh oc 
M° ®^ from Ferrarese (2002), which may be less fundamental than 
the spheroid-based relations (Kormendy & Ho 2013) but does hold 
empirically. For halo masses Mh < 10^^'^ M© at z = 0, the stel¬ 
lar mass scales with the halo mass as Mstar oc Mh'^, which using 
the Ferrarese (2002) relation implies Mstar oc M^'^. As discussed 
in §4.3, the presented model predicts Nqc oc Mstar = 3)°'^. 
Without a model that describes the further growth of galaxies and 
SMBHs between z ~ 3 and z = 0, there is no reason to assume 
that this relation can be directly connected to the observed Mstar- 
M* relation at z = 0. However, combining both expressions does 
yield Nqc oc M° in good agreement with the observed re¬ 
lation. In the context of this work, the good agreement between 
the A/^gc-M. relation in the z ~ 3 model and the z = 0 obser¬ 
vations suggests that the SMBH mass and the host galaxy stellar 
mass evolved roughly in proportion to each other between z ~ 3 
and z = 0 , as is indeed found in galaxy formation simulations (e.g. 
Sijacki et al. 2015). 

4.6 Summary of the presented results and comparison to 
previous models 

The results of this section show that the observed properties of GCs 
at z = 0 are reproduced by our model, indicating that GCs are 
fully consistent with being the products of the regular star and clus¬ 
ter formation process in the high-pressure environments of z > 2 
galaxies. These clusters were saved from tidal destruction in their 
dense natal environments by the redistribution of matter (and hence 
their ‘liberation’ into the galaxy halo) during hierarchical galaxy 
formation. The most important quantitative constraints from Fig¬ 
ures 3-5 are summarised in Table 2, which lists the GCMF peak 
mass Mpeak, the number of GCs Nqc, the specific frequency Tn, 
the ratio between the GC population mass and host galaxy stel¬ 
lar mass Mgg, tot/Mstar, and the ratio between the GC population 
mass and the host galaxy halo mass MGG,tot/Mh, at times tmerge 
and tnow (as defined in Figure 2). 

Previous models aimed at explaining the origin of GCs in the 
context of galaxy formation have almost all included a description 
for the origin of the GCMF (e.g. Prieto & Gnedin 2008; Katz & Ri- 
cotti 2014; Li & Gnedin 2014). These models reproduce the peaked 
shape just like the model presented here. However, they retrieve a 
universal peak mass by construction, because the dominant disrup¬ 
tion mechanism in these models is an environmentally-independent 
formulation of evaporation, in strong contrast with observations, 
theory, and simulations (e.g. Baumgardt & Makino 2003; Lamers, 
Gieles & Portegies Zwart 2005; Kruijssen et al. 2011). Our model 
does not make this assumption and instead reproduces the GCMF 
mainly by including the (dominant) evolutionary phase in the natal 
galaxy disc, while also accounting for the environmental depen¬ 
dence of the evaporation rate in the host galaxy halo. Owing to 
the important role of the rapid-disruption phase, the model does 
not predict a radial gradient of the peak mass or require extremely 
anisotropic GC orbits, both of which were predicted by earlier, 
evaporation-only models but are unobserved (e.g. Fall & Zhang 
2001; Vesperini et al. 2003). 

The dependence of the specific frequency on galaxy mass or 


© 2015 RAS, MNRAS 000, 1-29 


18 


/. M. D. Kruijssen 


Table 2. Properties of the GC populations at the times t = {tmerge, fnow}- 



Fornax (1) 

Fornax (2) 

[Fe/H] = -1.1 

[Fe/H] = -0.7 

[Fe/H] = -0.6 

Mpeak (Emerge) [10^ M©] 

1.2 

0.4 

2.7 

5.7 

7.8 

Mpeak(^now) [10^ M©] 

1.0 

0.3 

1.8 

2.3 

2.5 

-VGC (Emerge) 

1 

6 

23 

583 

902 
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metallicity is rarely considered in GC formation models. A recent 
example is presented in Katz & Ricotti (2014), who find a de¬ 
creasing number of GCs per unit host stellar mass. However, this 
is a direct consequence of their assumption of environmentally- 
independent GC disruption and their GC formation model, which 
assumes that GCs form in direct proportion to the dark matter halo 
mass and thus automatically represent a lower fraction of the stellar 
mass in galaxies with high stellar mass fractions. The specific fre¬ 
quency is therefore not predicted by the model, but follows directly 
from the assumptions. In model presented here, we do not make 
these assumptions, but predict the specific frequency using a phys¬ 
ical description in which GCs form from the ISM in galaxy discs 
{not in proportion to the halo mass) and undergo environmentally- 
dependent disruption, both of which are consistent with the variety 
of constraints presented in §2. 

Likewise, previous models reproduce the colour and metallic¬ 
ity distributions of GCs by assuming a relation between the number 
of GCs and the host galaxy stellar or halo mass, either in combina¬ 
tion with an environmentally-independent disruption model (e.g. 
Muratov & Gnedin 2010; Katz & Ricotti 2014; Li & Gnedin 2014) 
or without including a disruption model (Tonini 2013). In these 
models, the metal-poor GC population originates from the dwarf 
galaxies that are accreted during hierarchical galaxy growth, which 
have high specific frequencies and therefore lead to a peak at low 
metallicities, even though such a peak is not present in the field star 
metallicity distribution. Our model does not a priori assume any 
scaling between the number of GCs and the host galaxy proper¬ 
ties, and instead provides the physical explanation for these scaling 
relations. 

In view of the model’s simplicity, it is particularly compelling 
that it simultaneously reproduces a substantial number of observ¬ 
ables. This is encouraging and provides a fruitful starting point for 
future work, especially given the many predictions of the model, 
which are discussed in §5. 


5 PREDICTIONS AND IMPLICATIONS FOR 
GLOBULAR CLUSTER FORMATION AND 
EVOLUTION 

Having tested the two-phase model for GC formation against sev¬ 
eral of the key characteristics of GC populations in the local Uni¬ 
verse, it is appropriate to outline its implications, as well as to pro¬ 
vide quantitative predictions for future work. 


5.1 Most of the cluster disruption takes place during the 
rapid-disruption phase in the host galaxy disc 

As a function of both the initial and current GC masses. Figure 6 
shows the mass loss due to tidal shocking by GMCs during the 
rapid-disruption phase in the host galaxy disc relative to the to¬ 
tal dynamical mass loss integrated over the total GC lifetime. It is 
clear that the two-phase GC formation model predicts that most of 
the cluster disruption occurs in the host galaxy disc, i.e. during the 
rapid-disruption phase of Figure 2. During this phase, clusters re¬ 
side in the gas-rich component of the host galaxy and are disrupted 
by the tidal injection of energy during encounters with GMCs. All 
clusters with initial masses below Mi = 3 x 10^ M© are dis¬ 
rupted during this rapid-disruption phase, i.e. before they escape 
the host galaxy disc. For the most metal-rich clusters, the min¬ 
imum mass necessary to survive the rapid-disruption phase may 
have been as high as M ~ 10® M©. Only the presently most mas¬ 
sive (M Mpeak) or metal-rich ([Fe/H] ^ —0.7) GCs may have 
had such short dynamical times that they could respond adiabati- 
cally to the tidal perturbations, allowing them to experience most 
of their mass loss during the slow-disruption, evaporation-driven 
phase in the galaxy halo. 

5.2 Relations between the number of GCs and their host 
galaxies are in place at z = 1-2 

Because most of the cluster disruption takes place due to tidal 
shocking by GMCs in the host galaxy disc, we predict that the 
well-known relations between the number of GCs and their host 
galaxy properties (e.g. stellar mass, halo mass, or metallicity, see 
Figures 4 and 5) were already in place at z = 1-2, whereas the final 
masses and radii of GCs were attained by their quiescent evolution 
in galaxy haloes. Likewise, the GCMF at z = 1-2 should already 
have a peak mass similar to (or even higher than) the characteristic 
GC mass-scale seen at the present day, i.e. Mpeak ^ 10® ® M© for 
z < 2 (see Figure 3). This prediction strongly contrasts with mod¬ 
els explaining the peaked GCMF by relaxation-driven evaporation, 
for which the peak mass should decrease by 0.3 dex for every fac¬ 
tor of two in cosmic time (e.g. Fall & Zhang 2001), yielding peak 
masses of Mpeak = {10®, 10^ ®, 10^ ®} M© for z = {0.5,1, 2}, 
assuming a typical formation redshift of z = 3. The peak mass re¬ 
quired by the evaporation model at z = 2 is at least an order of 
magnitude lower than predicted by our model. 

By the time GCs have escaped into the host galaxy halo (which 
is typically after a few Gyr, see Table 1), they will have faded 
substantially compared to their initial luminosities. With HST, it 
is therefore not possible to detect the GC peak mass at z = 0.5- 
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Figure 6. This figure shows that the majority of cluster disruption takes place during the rapid-disruption phase in the host galaxy disc. Shown is the ratio 
between the dynamical mass loss during the rapid-disruption phase and the total dynamical mass loss AMcce/AMdyn,tot a function of the initial GC 
mass Mi (left panel) and the current GC mass M (right panel) for the five different galaxy models from Table 1 . The legend indicates the colour coding, 
where the {thin,thick} blue lines represent the {long,short}-fmerge Fornax model (towards the right and left in the figure, respectively). The grey-shaded area 
indicates where the dynamical mass loss during the rapid-disruption phase in the disc exceeds that during the slow-disruption phase in the halo. For reference, 
the upward-pointing arrows in the right-hand panel denote the peak masses Mpeak (which roughly coincide with the median masses) of the 2 ; = 0 GCMFs 
from Figure 3. Clusters with AMcce/AMdyn,tot = 1 in the left-hand panel are destroyed during the rapid-disruption phase, which occurs for the vast 
majority of the initial population. Omitting the short-fmerge Fornax model, only clusters with initial masses Mi ^ 10® M© (corresponding to current masses 
M Mpeak) or metal-rich, surviving GCs ([Fe/H] ^ —0.7) experienced most of their dynamical mass loss by evaporation in the host galaxy halo. 


2. Even with the James Webb Space Telescope (JWST), a signal- 
to-noise ratio of S/N > 1 may only be reached out io z — 0.5 
for massive (M > 10^ M©) GCs, requiring multi-hour expo¬ 
sures. Similar exposure times will be required for thirty-metre class 
telescopes such as the European Extremely Large Telescope (E- 
ELT) and the Thirty Meter Telescope (TMT). None the less, the 
intermediate-redshift evolution of the GCME is a critical discrimi¬ 
nator between the evaporation and two-phase GC formation models 
once observational facilities gain the required resolution and sensi¬ 
tivity. 


5.3 GCs were formed according to the same power-law 
ICMF as young clusters in the local Universe 

The prediction that the high, present-day GCME peak mass was al¬ 
ready reached dii z — 0.5-2 is not unique - the same prediction is 
made by models proposing that the ICME of GCs was peaked at 
birth (e.g. Parmentier & Gilmore 2007). However, in our model, it 
is assumed that the ICME follows the same power-law mass func¬ 
tion as seen for cluster formation in the local Universe. The good 
agreement of the modelled GC populations with those observed at 
z = 0 lends some support to this assumption - it suggests that GCs 
indeed formed according to the same power-law ICME as young 
clusters in the local Universe. As shown by Katz & Ricotti (2013), 
the YMC progenitors of GCs should be readily observable with 
JWST in the formation redshift range z = 2.3-8. The measured 
ICME should allow us to distinguish between the presented model, 
in which the ICME follows a power law, and those models propos¬ 
ing that the GC peak mass is an imprint of their formation process. 


5.4 The high-mass ends of the GMC and cluster mass 
functions at high redshift are set by the Toomre mass 

The two-phase GC formation model reproduces the properties of 
observed GC populations by assuming that the maximum mass- 
scale of the ICME is related to the maximum mass-scale for grav¬ 
itational instability (the Toomre mass) in the same way as holds in 
nearby star-forming galaxies. The associated prediction is that the 
high-mass end of the GMC mass function in z > 2 galaxies is set 
by the Toomre mass. With ALMA, it will be possible to observe the 
high-mass end of the cloud mass function in these galaxies, which 
will provide a direct test for the relation between the Toomre mass 
Mt and the GMC mass. In combination with the aforementioned 
measurement of the ICME (§5.3), it will also verify the adopted re¬ 
lation between the maximum GMC mass and the maximum YMC 
mass Me, which is assumed to satisfy Me = Mei,T = cFMt (see 
§3.1). 

5.5 The cosmic GC formation history peaks at a higher 
redshift than the cosmic SFH 

Galaxy mergers play a central role in our model, because they 
enable the migration of clusters out of the host galaxy disc into 
the galaxy halo. As has been discussed at length, this process 
strongly increases the long-term survival chances of the clusters. 
Because the merger rate increases with redshift, the cosmic GC 
formation history should peak at a higher redshift than the cos¬ 
mic star formation history. Earlier estimates by Peng et al. (2008, 
Eigure 20) suggest that the redshift difference is roughly /\z ~ 1 
(i.e. At ~ 1 Gyr for z = 2-4) and thus should be detectable given 
similarly-accurate GC age estimates. 
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5.6 There is a time delay between cluster formation and their 
migration into the galaxy halo 

Our model assumes a delay of several Gyr between cluster forma¬ 
tion in galaxy discs and their migration by mergers and/or stripping 
into the halo, which plays a critical role in determining the proper¬ 
ties of the GC population at z = 0. The good agreement between 
this simple model and the observations suggests that GC formation 
mechanisms that deposit GCs near-instantly in the halo are sub¬ 
dominant. As a result, either (1) the clusters that are formed during 
a merger itself have a low probability of escaping the host galaxy 
because they are destroyed on time-scales shorter than a dynami¬ 
cal time (see §2.2), or (2) the majority of all GC migration occurs 
through the tidal stripping of GCs from their host galaxies by more 
massive haloes, in which case migration is rarely associated with a 
starburst event. 

The first of the above points would suggest that all clusters that 
formed at least one galactic dynamical time before a major merger 
can safely migrate into the halo. Observations of the galaxy merger 
pair NGC 4038/39 (the Antennae galaxies) by Bastian et al. (2009) 
show that surviving clusters without disc kinematics (indicative of 
migration) indeed have ages of 200-500 Myr, which is confirmed 
by numerical simulations of galaxy mergers (Kruijssen 2011). This 
indicates a formation time that precedes the most recent encounter 
by > 100 Myr and thus supports point (1). The second point can 
be tested using the mass ratio-dependence of the merger rate in the 
Millennium Simulation. Genel et al. (2009) show that at z ~ 3, 
only galaxies with halo masses Mh < 10^^ M© typically merge 
with more massive galaxies and thus lose their GCs by stripping. 
At higher halo masses (corresponding to [Fe/H] > —0.8, see Ta¬ 
ble 1), the majority of the merging takes place with haloes of similar 
masses, indicating that the migration of GCs by the tidal stripping 
of their host galaxies is not significant. This metallicity cut range 
covers a minority of the GC metallicity range, confirming that point 
(2) above typically holds too. 


5.7 The increase of the specific frequency at high galaxy 
masses is caused by a combination of satellite galaxy 
accretion and merger-induced cluster formation 

We see that GCs with metallicities [Fe/H] > —0.8 have migrated 
into the halo by major mergers, which is possibly associated with 
a starburst if the host galaxy is gas-rich, whereas lower-metallicity 
GCs have migrated by tidal stripping. In the Virgo Cluster, the mean 
metallicities of the GC populations reach [Fe/H] ~ —0.8 for a 
host galaxy magnitude Mb ~ —20 (or a stellar mass Mstar ~ 
10^^ M©), coinciding exactly with the magnitude at which the spe¬ 
cific frequency Tn reaches a minimum (see Figure 6 of Peng et al. 
2008). As discussed in §4.3, it is likely that the increase of Tn in 
more massive galaxies (Mstar >10^^ M©) at least partially arises 
due to the accretion of GCs from lower-mass galaxies, but an ad¬ 
ditional contribution may come the change of GC migration agent 
discussed here. If the migration of clusters in haloes with masses 
Mh > 10^^ M© is associated with a starburst in which large num¬ 
bers of clusters are born, at least some of the newly-formed clusters 
could find safety in numbers and survive despite the low survival 
probability of clusters born at the peak SFR of galaxy mergers. 

Future work can quantify the relative contributions of both 
mechanisms for increasing Tat at high galaxy masses with a num¬ 
ber of simple measurements. Firstly, the mass fraction accreted in 
the form of dwarf galaxies (which have high Tn) should increase 
with galaxy mass more steeply than Mace,tot oc M^iir to offset the 


intrinsic decrease of Tn with galaxy mass (see Figure 4). Secondly, 
the metallicities of the GCs leading to the increase of the specific 
frequency at high galaxy masses differs strongly for both mecha¬ 
nisms. If tidal stripping dominates, the extra GCs should predom¬ 
inantly be metal-poor and the increase of the mean GC metallicity 
with galaxy mass should diminish. If the migration event is asso¬ 
ciated with a starburst that supplies additional GCs, the extra GCs 
should be metal-rich and the correlation between mean GC metal¬ 
licity and galaxy mass should steepen. Figure 9 of Peng et al. (2008) 
shows that Tn increases both for the metal-rich and the metal-poor 
GCs for galaxies with Mb < —20, suggesting that both mecha¬ 
nisms should be at work. 


5.8 The specific frequency increases with galaxy clustering 

The galaxy merger time-scale varies with the cosmic environment, 
e.g. between highly clustered and (near-)isolated field galaxies. A 
general prediction of the presented model is that the specific fre¬ 
quencies of field galaxies (which have low merger rates) should be 
systematically lower than those of galaxies in dense environments 
(which have high merger rates). This environmental dependence of 
Tn arises because the modelled duration of the rapid-disruption 
phase is set by the galaxy merger time-scale. In other words, the 
specific frequency of otherwise identical progenitor galaxies with 
different merger time-scales should scale as Tat oc tm^ge- Obser- 
vationally, this variation can be constrained by considering GCs 
at fixed metallicity, which removes the dependence of Tn on the 
galaxy mass at GC formation seen in Figure 4. Indeed, Peng et al. 
(2008) find an increase of the specific frequency towards the centre 
of the Virgo Cluster at fixed g — z colour (i.e. at fixed GC metallic¬ 
ity), which at least qualitatively confirms the dependence of the GC 
survival fraction on the galaxy merger rate that is predicted by the 
proposed, two-phase model. Note that even if galaxy mergers are 
not responsible for the migration of GCs into a gas-poor environ¬ 
ment, but instead the gas is removed (e.g. by feedback, strangula¬ 
tion, or ram pressure stripping, van den Bosch et al. 2008), then the 
enhanced quenching of star formation activity observed in galaxy 
clusters (‘environment quenching’, see e.g. Peng et al. 2010) sug¬ 
gests that the rapid-disruption phase may have a shorter duration in 
clustered environments, (although the time-scales may be of a sim¬ 
ilar order of magnitude, for e.g. strangulation see Peng, Maiolino 
& Cochrane 2015). In addition, such early quenching of star for¬ 
mation prevents the further growth of the field star population after 
the majority of the GCs has formed, which further boosts the spe¬ 
cific frequency. All of these mechanisms should lead to a trend of 
an increasing specific frequency with galaxy clustering. 


5.9 The GC evaporation time-scale decreases with the GC 
metallicity 

To obtain a near-universal GCMF, our model requires a weak nega¬ 
tive metallicity gradient of a galaxy’s GC population, such that the 
disruption time-scale due to evaporation decreases with the metal¬ 
licity (see Table 1 and Figure 3). In §4.2, we showed that the re¬ 
quired gradient is quantitatively consistent with that observed for 
GCs in the Milky Way. It was proposed that this gradient is an im¬ 
print of the galaxy mass (and hence metallicity) dependent redistri¬ 
bution of GCs in mergers and/or tidal stripping, in which metal-rich 
GCs formed at higher binding energies (and experience stronger 
tidal fields). If true, then the relation between [Fe/H] and the dis¬ 
ruption time-scale due to evaporation provided in equations (15) 
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and (16) may hold universally in galaxies where dynamical friction 
on GCs is negligible (also see §5.12). 


5.10 Fossil stellar streams in the halo from evaporated GCs 
are more metal-rich than the surviving GCs 

Observational missions and surveys such as Gaia, SEGUE and 
RAVE are capable of identifying and characterising the remnants of 
tidally disrupted stellar systems in the Galactic halo, be it tidal de¬ 
bris from disrupted dwarf galaxies or from evaporated GCs. It was 
shown in §5.1 and the right-hand panel of Eigure 6 that for metal- 
licities [Fe/H] < — 0.7 most of the mass loss takes place during the 
rapid-disruption phase in the host galaxy disc, implying that little 
coherent debris of these clusters should be preserved in the Galac¬ 
tic halo. By contrast, metal-rich GCs ([Fe/H] ^ —0.7) should typi¬ 
cally have lost most of their mass by evaporation in the halo. This is 
a corollary of §5.9, which states that the GC evaporation time-scale 
decreases with metallicity. As we will discuss in §5.13, GCs with 
metallicities [Fe/H] ^ —0.7 and typical masses (M ~ 10^ M©) 
may have lost 90% of their initial mass, in strong contrast with the 
substantially lower mass loss (^ 50%) of lower-metallicity GCs 
at the same present-day mass. As a result, the tidal debris of dis¬ 
rupted GCs in galaxy haloes should be dominated by metallicities 
larger than the median GC metallicity. With Gaia and its spectro¬ 
scopic companion surveys (e.g. Gaia-ESO), it will be possible to 
systematically test this prediction for the Galactic halo. 


5.11 The radii of GCs are set by structural evolution and are 
unrelated to GC formation or disruption 

While our model predicts that most of the cluster disruption oc¬ 
curs during the disc phase of Eigure 2, most of the GCs’ structural 
evolution takes place during the halo phase. As a result, the empir¬ 
ical relation seen in present-day GC populations between the GC 
half-mass density and GC (peak) mass (M oc see e.g. Eall & 
Rees 1977; Gnedin & Ostriker 1997; McLaughlin & Eall 2008) is 
caused by the evolution of (primarily) the GC radius (and to a lim¬ 
ited degree the GC mass) in the host galaxy halo. The variation of 
the GC mass with the half-mass density is therefore unrelated to the 
conditions at GC formation. This is to be expected theoretically - 
unless initially rn oc with (3 < —1/3 (at odds with theory and 
observations, see §2.1.1), the half-mass relaxation time increases 
with the cluster mass, causing low-mass clusters to expand more 
rapidly than massive clusters. After a time thaio of free expansion 
in the host galaxy halo, this leads to a relation between the half- 
mass density and the cluster mass of M oc p^l (Gieles, Heggie 
& Zhao 2011), which is roughly consistent with the increase of 
the peak mass with the half-mass density from McLaughlin & Eall 
(2008). The observation that the typical GC mass increases with the 
half-mass density is thus not causal in the sense that high-density 
GCs dissolve more rapidly, but simply arises due to their structural 
evolution (Gieles, Heggie & Zhao 2011). 


5.12 Dynamical friction has affected only a minority of the 
GC population 

Given the properties of the host galaxy, a minimum cluster mass- 
scale can be derived above which the dynamical friction time-scale 
to the centre of the galaxy is shorter than the host galaxy merger 
time-scale, in which case the clusters are destroyed by dynamical 


friction before they can escape into the halo. Using equation (6) of 
Kruijssen (2014), we obtain 


Mdf,min = 10® M© 


( Emerge 

V2^ 



_f_) 

200 km s ^ / 


.(18) 


Setting the circular velocity equal to the virial velocity of the halo 
at z = 3, we have V — , where H{z) is the 

Hubble constant at redshift z, i.e. H{z) — iTo[^m(l + 2 ;)^ + (l — 
^m —^a)(1+^)^+^a]^^^, which in combination with the angular 
velocity U as a function of stellar (and hence halo) mass from equa¬ 
tion (6) also defines the radius R — U/U as a function of galaxy 
mass. Substitution into equation (18) then yields minimum mass- 
scales for dynamical friction into the galaxy centre® ranging from 
10®-® M© for the late-merger Eornax model to 10® ®-10® ® M© 
for all other galaxy models, indicating that dynamical friction may 
have affected a small part of the GC population (thereby possibly 
fueling the formation of nuclear clusters), but never dominated dur¬ 
ing early GC evolution. This confirms the result of Jordan et al. 
(2007), who found that dynamical friction cannot explain the expo¬ 
nential truncation mass at the high-mass end of the GCME. 


5.13 Typical GCs were only 2-3 times more massive at birth 

Considering that the two-phase GC formation model reproduces 
several of the observed properties of GC populations, it may be 
used to constrain GC formation theories aiming to reproduce the 
chemical abundance patterns seen in Galactic GCs. In §2.1.3, we 
very briefiy summarised some of the main assumptions, predic¬ 
tions, and problems of these models. One of these is the mass bud¬ 
get problem, requiring GCs to have been much more massive at 
birth than they are today, i.e. by a factor of 10-100 (for explaining 
the Na-0 anti-correlation in the AGB enrichment model) or by a 
factor of > 20 (for forming a sufficiently massive enriched stel¬ 
lar population in the ERMS enrichment model). The two-phase GC 
formation model allows us to predict initial-to-current mass ratios 
Mi/M for GCs across the observed GC mass and metallicity range 
and test the assumptions of these enrichment models. 

Eigure 7 shows Mi/M as a function of the present-day GC 
mass M for the five galaxy models of Table 1 at the time of GC 
migration and at the present day (tmerge and tnow as defined in 
Eigure 2). We see that the initial-to-current mass ratio is a function 
of the current GC properties, decreasing strongly with the GC mass 
and increasing with the metallicity. There is little evolution between 
Emerge and tnow, which again indicates that most of the disruption 
took place in the GCs’ high-redshift birth environment. Only low- 
mass and/or metal-rich GCs may have lost > 90% of their initial 
mass. If we assume the simple case that both the AGB and ERMS 
enrichment scenarios require GCs to have been at least 10 times 
more massive at birth (which is a strict lower limit for both models, 
see §2.1.3), then Eigure 7 provides an upper limit to the present-day 
GC mass Mxio,max below which the AGB or ERMS enrichment 
models may have operated. This mass limit is typically much lower 
than the peak mass of the GCME. Only at metallicities [Fe/H] ^ 
-0.7 we find MxlO,max ~ Mpeak- 

At metallicities [Fe/H] < —0.7 the majority of the GC popu¬ 
lation has a mass loss history that is inconsistent with the AGB or 
ERMS enrichment models, as most GCs were much less than 10 


® To calculate these numbers, we use the Planck 2013 cosmological param¬ 
eters, with Ho = 67.3 km Mpc“^, Um = 0.315 and Ua = 0.685 
(Planck Collaboration et al. 2014). 
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Figure 7. This figure shows that only a minority of the present-day GC pop¬ 
ulation satisfies the requirement of the AGB and FRMS enrichment models 
that GCs were > 10 times more massive at birth. Shown is the initial-to- 
current GC mass ratio Mi/M as a function of the current GC mass M 
for the five different galaxies from Table 1 . The legend indicates the colour 
coding, where the {thin,thick} blue lines represent the {long,short}-f merge 
Fornax model (towards the top right and bottom left in the figure, respec¬ 
tively). Dashed lines show the mass ratio at the end of the rapid-disruption 
phase (fmerge) and solid lines show the mass ratio at 2 ; = 0 (tnow)- The 
horizontal dotted line represents zero mass loss (Mi/M = 1) and the grey- 
shaded area indicates where the initial mass is > 10 times larger than the 
current GC mass. For comparison to the solid lines, the upward-pointing 
arrows denote the peak masses Mpeak (which roughly coincide with the 
median masses) of the 2 ; = 0 GCMFs from Figure 3. The condition that 
GCs were >10 times more massive at birth is only satisfied up to a max¬ 
imum mass Mxio,max, which increases from Mxio,max ~ 10^-^ Mq 
for the short-fmerge Fornax model to Mxio,max ~ 10^-^ Mq for the 
[Fe/H] = —0.6 model. Note that Mxio,max < ^peak for roost of the 
metallicity range ([Fe/H] < —0.6), showing that the majority of GCs had 
initial masses only a few times larger than their present-day masses. 

times more massive at birth. As explained in §2.1.3, both enrich¬ 
ment models assume that most of the unenriched stars are lost in 
order to avoid the mass budget problem. In other words, they pre¬ 
dict that the enriched fraction must increase with the initial-to-final 
mass ratio, because it reflects the amount of mass loss. Combining 
this prediction with the results of Figure 7, we see that according 
to the AGB and FRMS models the enrichment fraction should de¬ 
crease with the GC mass and increase with the metallicity, becom¬ 
ing very small for M > 10^ M© and [Fe/H] < —0.7. The contrary 
is observed - as discussed in §2.1.3, chemically enriched stellar 
populations have been observed in all considered GCs with masses 
10^ ^ M/M© ^ 10®-^ and metallicities -2.3 < [Fe/H] < -0.3, 
and the enriched fraction increases with the GC mass (see e.g. Fig¬ 
ure 2 of Conroy 2012). 

It was shown in §4 that the predicted initial-to-current mass ra¬ 
tios of Figure 7 are consistent with the observed ratios between the 
number of GCs and the field star population (i.e. the specific fre¬ 
quency). If the GCs would have attained their present-day masses 
by losing more mass, then a larger field star population would have 
been necessary. In other words, the host galaxies simply do not 
contain enough stellar mass for the GCs to have lost more mass 
than predicted by the two-phase model, even though at first sight 
the small fraction of the total stellar halo mass constituted by GCs 
may suggest otherwise (see §5.14). As we have seen in §2 and 
§3, the reason for the low GC mass fraction in galaxy haloes is 


that during a star formation event, only some fraction of the stel¬ 
lar mass ends up in the clusters with masses M > 10 ® M 0 that 
could survive to become GCs. As a result, galaxy haloes through¬ 
out the Universe would need to contain several times more field 
stars than observed to be consistent with the AGB and FRMS en¬ 
richment models, unless the stellar initial mass function, the CFE 
or the ICMF were altered, in conflict with theory and observations. 
It is therefore more likely that the fraction of enriched stars reflects 
some aspect of the conditions at GC formation (as it does in the 
EDA enrichment model) than that these stars formed from AGB or 
FRMS ejecta in a second star formation episode. 


5.14 The minimum mass for the formation of multiple stellar 
populations in GCs was ~ 10^ M© sit z > 2 

All GC enrichment models discussed in §2.1.3 suggest that the 
presence of a chemically enriched stellar population requires some 
minimum cluster mass (or density, which for oc with 
/3 < 1/3 can still be translated to a minimum mass). By comparing 
the total enriched mass in the Galactic stellar halo to that in GCs, it 
is possible to empirically constrain this minimum mass for cluster 
enrichment. About 2% of the field star population in the Galac¬ 
tic halo consists of enriched stars (e.g. Carretta et al. 2010, also see 
Vesperini et al. 2010). This corresponds to a mass of ~ 2 x 10^ M© 
(e.g. Bell et al. 2008), similar to the total enriched stellar mass 
presently in GCs (also 1-2 x 10^ M©, assuming that presently 
50% of the GC stars are enriched). As such, a total of ~ 4% of 
all Galactic halo stars are enriched. If we assume that enrichment 
only took place in clusters above some minimum initial mass, we 
can use the mass loss histories provided by our two-phase GC for¬ 
mation model to determine the minimum mass for GC enrichment 
to produce a stellar halo in which a total of 4% of all stars are en¬ 
riched.® 

To address the origin of the enriched halo stellar population, 
we first determine what fraction /gc,! of the entire coeval stellar 
population is constituted by stars that were born in (but possibly 
lost from) the present, surviving GCs. For the five galaxy models 
considered here, this number varies in the range /GC,i = 2-16%, 
with a median of /Gc,i = 7% - the other 93% was never part of 
the GCs that survived to the present day. Given that the present GC 
enrichment ratio (enriched : unenriched ~ 1 : 1 ) is representa¬ 
tive of the initial enrichment fraction in GCs to within a factor of 
2 (Decressin, Baumgardt & Kroupa 2008), a fraction 1/3-2/3 of 
the lost stars may be enriched, i.e. 1/3-2/3 x 7% ~ 2-5% of the 
stellar halo consists of enriched stars that once resided in surviving 
GCs. This matches the observed enriched fraction in the Galactic 
halo and therefore requires that all of the enriched stars have been 


® The answer to this question necessarily relies on the relative escape rates 
of stars from the enriched and non-enriched populations. If non-enriched 
and enriched stars have been lost from GCs at a similar rate throughout cos¬ 
mic history, then the similar total enriched stellar masses in the halo and in 
GCs suggests that enriched stars formed in GCs above some minimum mass 
such that on average they lost half of their mass through dynamical mass 
loss. In principle, this minimum mass could be estimated for each of the 
model galaxies considered here. However, the three GC enrichment models 
also predict that the enriched population should be more centrally concen¬ 
trated than the non-enriched population. This chemical segregation causes 
non-enriched stars to be lost preferentially, which may increase the ratio 
of enriched to non-enriched stars by up to a factor of two (e.g. Decressin, 
Baumgardt & Kroupa 2008). 
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supplied by the present GCs. The minimum initial mass of the sur¬ 
viving GCs Mi,min should thus reflect the minimum cluster mass 
required for enrichment Menr,min. 

Our five galaxy models yield Mi,min = Menr,min = 
Mq (see Figure 6) with a median of Menr,min = 
10^'^ Mq/° This high mass is required by the small fraction of 
enriched stars in the Galactic halo, and is uncertain by a factor of 
~ 2 depending on how representative the present GC enrichment 
fraction is of the initial GC population (Decressin, Baumgardt & 
Kroupa 2008). Considering that most of the surviving GCs had ini¬ 
tial masses exceeding this limit, it is possible that multiple stellar 
populations are indeed a universal property of GCs. 


5.15 Most GCs could not have survived the high gas densities 
required for pristine gas accretion 

Finally, we turn to the gas accretion problem of the AGB enrich¬ 
ment model, which requires a pristine gas mass comparable to the 
initial GC mass to be accreted for the formation of the enriched 
stellar population (see §2.1.3). A possible mechanism for supply¬ 
ing the gas (e.g. Conroy & Spergel 2011; Maxwell et al. 2014) is 
Bondi-Hoyle-Lyttleton (BHL) accretion (Hoyle & Lyttleton 1941; 
Bondi & Hoyle 1944; Bondi 1952), which occurs at a rate 

(dM\ _ 2aG'^M^pisM 

where a is a constant of order unity, Urei is the relative velocity 
between the accreting body and the surrounding material, and Cg is 
the sound speed. One can define a time-scale for accreting the ini¬ 
tial mass of the accretor by writing tsHL = M/(dM/dt)BHL oc 
(where we have assumed Vrei ^ Cg). This expres¬ 
sion shows that accretion proceeds more rapidly in a higher-density 
ISM. However, the disruption time-scale due to tidal shocks has a 
very similar scaling of tcce oc Mpjg^^, indicating that a higher- 
density ISM disrupts clusters more efficiently. Because the density- 
dependence of the disruption time-scale is steeper than that of the 
accretion time-scale, there must exist a maximum ISM density be¬ 
low which significant gas accretion can occur before a cluster of 
mass M is disrupted, i.e. tsHL/tcce oc < 1. Vice versa, 

given an ISM density there exists a minimum mass above which a 
cluster can survive long enough to accrete. 

The competition between accretion and disruption is visu¬ 
alised in Figure 8, which shows the time-scale ratio tsHL/Ace 
as a function of the initial (dotted) and present-day (solid) GC 
mass for the five galaxy models of Table 1. For the relative ve¬ 
locity we have adopted the circular velocity V (rather than a 
velocity dispersion), which is appropriate for clusters on plung¬ 
ing, eccentric orbits (as proposed by Maxwell et al. 2014). Fig¬ 
ure 8 shows that the conditions necessary for significant BHL 
accretion are so disruptive that only clusters with very high ini¬ 
tial masses may have survived the process. Pristine gas accretion 


If enrichment occurs in clusters above some critical half-mass density 
rather than above some critical mass, it is possible that the corresponding 
mass limit is higher in the local Universe than it was at the time of GC 
formation. The initial radii of YMCs may increase with decreasing pressure 
or ambient density (see §2), and hence with cosmic time. As explained in 
§3, we have scaled the YMC radius by the mean tidal density, leading us 
to adopt radii at 2 ; = 3 a factor of 2-3 smaller than those seen at 2 ; = 0. 
In that specific case, the present-day threshold mass for the formation of an 
enriched population would currently be ~ 10 times larger than at 2 : = 3, 
i.e. ^ 10® M©. 



M [Mo] 

Figure 8. This figure shows that only a minority of the present-day GC pop¬ 
ulation can satisfy the requirement of the AGB enrichment model that GCs 
should accrete a pristine gas mass similar to their initial mass. Shown is the 
time-scale ratio between the Bondi-Hoyle-Lyttleton (BHL) accretion time- 
scale and the disruption time-scale during the disc phase (see Figure 2), as 
a function of the initial (dotted lines) and current (solid lines) GC masses 
for the five different galaxies from Table 1 . The legend indicates the colour 
coding, where the {thin,thick} blue lines represent the {long,short}-tmerge 
Fornax model (towards the bottom and top in the figure, respectively). The 
grey-shaded area indicates where ^ehl < tcce, which is required for gas 
accretion to outpace cluster disruption by the high-density environment. For 
comparison to the solid lines, the upward-pointing arrows denote the peak 
masses Mpeak (which roughly coincide with the median masses) of the 
2 ; = 0 GCMFs from Figure 3. The condition for pristine gas accretion 
is only satisfied above some minimum mass Mace,min that exceeds the 
GCMF peak mass (Mace,min > Mpeak) for almost all metallicities. 


without disrupting the cluster is not possible below a mass-scale 
Mace,min, which incrcascs with metallicity from {initial, current} 
masses of Mace,min ~ {10^'^, 10^'®} Mq at [Fe/H] = —2.3 to 
Maee,min ~ {10® ®, 10®'^} M© at [Fo/H] = —0.6. Bccausc this 
analysis strictly concerns the dynamics of gas accretion and tidal 
disruption while excluding any form of feedback, these masses rep¬ 
resent firm lower limits for pristine gas accretion. 

Comparing Maee,min(tnow) to Mpeak across the galaxy mass 
range shows that Maee,min > Mpeak for all galaxies other than 
the long-tmerge Fomax model, indicating that a minority of young 
GCs can have accreted pristine gas. Moreover, Mace,min > Me 
in all galaxies other than the Fornax models, suggesting that few 
clusters massive enough to undergo significant pristine gas accre¬ 
tion ever formed at all. We conclude that high-pressure environ¬ 
ments are generally unable to drive the accretion of pristine gas 
onto young GCs without disrupting them. 

At the extremely high GC masses where BHL accretion could 
possibly proceed without disrupting the cluster, the results of Fig¬ 
ure 8 imply that the amount of dilution by pristine gas in the AGB 
enrichment model should increase with the GC mass. In the AGB 
enrichment model, it is this dilution which turns the Na-0 corre¬ 
lation predicted by stellar evolution models into the observed anti¬ 
correlation. As we have seen, metal-rich galaxies at z > 2 are more 
disruptive to their young GCs than metal-poor ones, hence Figure 8 
shows that the Na-0 anti-correlation should weaken towards in¬ 
creasing metallicity, because the minimum mass for dilution should 
increase. This suggests that metal-rich GCs like NGC 104, 6388, 
and 6441 should only have been able to develop a weak Na-0 anti- 
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correlation (or none at all) if the AGB enrichment model holds. 
However, no such dependence on metallicity is observed (e.g. Grat- 
ton, Carretta & Bragaglia 2012). 

Finally, comparing Mace,min to the maximum mass 
Mxio,max below which GCs were > 10 times more massive 
at birth (see Figure 7), it is clear that Mace,min /Mxio, max ^ 4 
for all galaxies, indicating that the two conditions required by the 
AGB enrichment model (allowing substantial mass loss as well 
as pristine gas accretion) cannot be satisfied simultaneously. This 
presents a serious challenge to AGB enrichment as the explanation 
for the observed abundance anomalies. 


5.16 Summary of predicted critical GC mass scales 

Throughout this section, we have discussed several predictions of 
the two-phase GC formation model for future observational tests, 
as well as its implications for existing GC formation and evolu¬ 
tion models. Among these were several minimum and maximum 
mass-scales above or below which several physical mechanisms 
could operate. For reference, these mass limits are summarised in 
Table 3, which shows the minimum initial GC mass of the present- 
day GCs Mi,min, the minimum GC mass for dynamical friction to 
the galactic centre during the disc phase Mdf,min, the minimum 
initial (tform) and current (tnow) GC masses needed for the accre¬ 
tion of pristine gas during the disc phase Mace,min (as required 
by the AGB enrichment model), and the maximum current GC 
mass below which the initial mass was >10 times higher at birth 
Mxio,max (as required by the AGB and FRMS enrichment mod¬ 
els). Note that for all models Macc,min > Mxio,max, indicating 
that no GC can both have accreted a substantial amount of pristine 
gas and have been >10 times more massive at birth. 


6 OPEN QUESTIONS AND FUTURE DIRECTIONS OF 
RESEARCH 

In §3, we used the observational constraints on GC formation from 
§2.1 and the theoretical state of the field from §2.2 to formulate a 
two-phase model for GC formation. While this model reproduces 
the main observed properties of GC populations and makes several 
predictions for future work, its simplicity should not be underesti¬ 
mated. The model describes the big picture of GC formation and 
evolution until the present day, but there are several open questions 
that still require answering. 

(i) Based on the pronounced relation between the GMC mass 
function and the ICMF of star-forming galaxies in the local Uni¬ 
verse (see §2.1.1), we have connected the maximum cluster mass 
to the Toomre mass, i.e. the maximum mass-scale below which 
self-gravity can overcome galactic shear. While this is a reasonable 
step to make and has been tested in the nearby Universe (Kruijssen 
2014; Adamo et al. 2015), observational confirmation is still lack¬ 
ing at high redshift. Are the Toomre mass and the maximum cluster 
mass related in the same way at z > 2 as they are at z = 0? If so, 
what are the masses of the most massive GCs that ever formed? 
The most massive clumps seen ‘m z > 2 galaxies have masses of 
~ 10® M©, consistent with the Toomre mass in these systems (e.g. 
Genzel et al. 2011; Swinbank et al. 2012). While these clumps may 
represent the extreme end of the mass-scale and may not be the 
most abundant, the obvious question is whether they ever formed 
any YMCs with masses M > 10^ M© (i.e. more massive than the 


initial masses of most GCs), and whether any such YMCs then es¬ 
caped the gas-rich component of their host galaxy. Above some 
mass-scale, the clumps and YMCs may have spiralled in to the 
galaxy’s centre by dynamical friction, but we have seen in §5 that 
the required minimum mass is well in excess of the exponential 
truncation mass of the ICMF. None the less, could dynamical fric¬ 
tion have moved extreme YMCs into different environments? In the 
coming years, these questions should be addressable at z > 2 us¬ 
ing ALMA to study massive gas clumps and large-scale telescopes 
equipped with adaptive optics (e.g. E-ELT, TMT) to study YMCs. 

(ii) A key ingredient for reproducing the observed properties of 
GC populations in nearby galaxies is an increasing ambient ISM 
density pisM (or more fundamentally its pressure) with host galaxy 
mass Mstar during the formation of the GC progenitor clusters. 
This is required for obtaining the observed dependence of the spe¬ 
cific frequency on the metallicity and galaxy stellar mass at for¬ 
mation - this relation between pisM and Mstar sets the disruption 
time-scale during the rapid-disruption phase tece and thus deter¬ 
mines what fraction of the initial cluster population survives until 
they migrate into the halo by a galaxy merger and/or tidal stripping. 
How universal is the pisM-Mstar relation? Again, using ALMA to 
map the molecular gas properties of z > 2 galaxies will be instru¬ 
mental in answering this question. 

(iii) In the presented model, galaxy mergers and the tidal strip¬ 
ping of accreted satellites facilitate the migration of GC-like YMCs 
from the gas-rich host galaxy disc into the galaxy halo. In the con¬ 
text of hierarchical galaxy formation, it seems indisputable that this 
process takes place. However, it is very well possible that addi¬ 
tional migration mechanisms exist. Eor instance, secular evolution¬ 
ary processes such as dynamical interactions with the gas clumps 
seen in high-redshift galaxies could eject massive clusters from 
the disc. If the Galactic disc was in place before the formation 
of the bulge (e.g. Bournaud, Elmegreen & Elmegreen 2007), this 
scenario could explain why the metal-rich GCs in the Milky Way 
exhibit a net co-rotation with the disc (Erenk & White 1980). If 
multiple processes contribute to the ‘liberation’ of GCs into their 
host galaxy haloes, the key question is what the relative contribu¬ 
tions and characteristic time-scales of these processes are. In the 
absence of clumpy galaxies in the local Universe to provide any 
observational constraints, the most promising avenue to address 
this question will be to quantify the migration mechanisms in high- 
resolution galaxy formation simulations with sufficient resolution 
to resolve the formation of gravitationally bound, massive clusters. 

(iv) As discussed in §2.2, it is well-understood how GC radii 
evolve in the quiescent environment of galaxy haloes. However, it 
is unclear how the structural evolution of YMCs proceeds in the 
erratic tidal fields of the high-pressure galaxy discs in which they 
formed. During the rapid-disruption phase of the two-phase model 
(see Eigure 2), the initial cluster radii and their time-evolution influ¬ 
ence how susceptible the clusters are to disruption by tidal shocks 
(for which the disruption time-scale increases with the cluster’s 
volume density). At the same time, tidal shocks themselves may 
induce cluster compression or expansion, suggesting that cluster 
disruption could be regulated (if tidal shocks cause compression) 
or could be a runaway process (if tidal shocks cause expansion). 
In either case, cluster survival is intimately linked not only to their 
mass evolution, but also to their radius evolution. A related question 
is whether there exists a maximum stellar volume density above 
which two-body relaxation acts as a spring that drives expansion 
and stops further compression. If so, this could lead to a depen¬ 
dence of the cluster radius on the cluster mass of rn oc 
(resulting in constant densities) for masses above which the maxi- 
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Table 3. Predicted GC mass limits for surviving GCs, dynamical friction, gas accretion, and > 90% mass loss. 



Fornax (1) 

Fornax (2) 

[Fe/H] = -1.1 

[Fe/H] = -0.7 

[Fe/H] = -0.6 

[105 M©] 

1.32 

0.29 

4.08 

10.73 

17.24 

-^df,min(fform) [10^ M©] 

6.82 

94.88 

44.24 

332.79 

1938.50 

-^acc,min(fform) [10^ M©] 

1.68 

1.68 

6.24 

19.60 

41.44 

-^acc,min(fnow) [10^ M©] 

0.64 

0.97 

2.63 

5.15 

12.42 

-^X 10,max(fnow) [10^ M©] 

0.14 

0.03 

0.43 

1.31 

2.27 


mum volume density is reached. Depending on what the threshold 
density would be, this could limit the survival of massive clusters 
during the rapid-disruption phase since tcce oc ph- 

(v) The simple model presented here is fundamentally unable to 
address the spatial evolution of cluster populations beyond the sim¬ 
ple classification of clusters residing in the host galaxy disc or halo. 
How would GC populations evolve spatially under the conditions 
described in the proposed model? While this question should ide¬ 
ally be addressed using numerical simulations, it is possible to get a 
rough idea by considering the current literature. There exists a wide 
range of models in which the GCs are initialised in the galaxy halo 
(i.e. at tmerge in our model, see Figure 2) and lose mass (largely) 
independently of the galactic environment (e.g. Chernoff & Wein¬ 
berg 1990; McLaughlin & Fall 2008; Prieto & Gnedin 2008; Li 
& Gnedin 2014). Such models can reproduce the observed GCMF, 
colour bimodality, and the spatial distribution of GCs (e.g. Prieto & 
Gnedin 2008; Li & Gnedin 2014). Although the cluster disruption 
model used here is environmentally dependent in general, it is most 
strongly so during the rapid-disruption phase in the host galaxy 
disc. Given that very few GCs are destroyed during the halo phase, 
they should be (re)distributed during hierarchical galaxy growth in 
the same way as in environmentally-independent models. This sug¬ 
gests that the two-phase GC formation model should also yield spa¬ 
tial distributions at z = 0 that are consistent with those observed. 

(vi) Thus far, we have considered massive cluster formation at a 
particular redshift z = 3 taken to be representative of the formation 
redshift of present-day GCs. How do the model predictions change 
when accounting for a realistic spread of GC formation redshifts 
and progenitor galaxy metallicities? Or conversely, what is the ex¬ 
act redshift range over which most of the GCs formed in our model? 
Some models have required GC formation to take place at z = 7- 
15, i.e. near the epoch of reionization, based on their present-day 
spatial distribution and hence their (initial) binding energies (e.g. 
Moore et al. 2006; Brodie & Strader 2006). However, these models 
did not account for the two-phase formation and disruption physics 
presented here and therefore they have limited information on the 
binding energies at which (the surviving) GCs actually formed. In¬ 
cluding the appropriate physics should change the inferred range of 
GC formation redshifts. 

(vii) A related question is how the formation redshift range 
varies with the host galaxy mass or GC metallicity. Many mod¬ 
els assume that metal-poor GCs formed in low-mass galaxies be¬ 
fore the metal-rich ones did in massive galaxies (e.g. Muratov & 
Gnedin 2010). A similar relation is suggested by relative age mea¬ 
surements of GCs (e.g. Forbes & Bridges 2010). However, the star 
formation histories of low-mass galaxies peak later in cosmic time 
than those of massive galaxies (often referred to as ‘downsizing’, 
e.g. Cowie et al. 1996). How can these two observations be united? 
Did metal-poor GCs indeed form at higher redshift than metal-rich 
GCs? Or do they trace the peak star formation activity of their host 
galaxies? Because the merger (and hence GC ‘liberation’) rate in¬ 


creases with redshift, the peak of the GC formation history should 
be shifted towards higher redshifts than the peak star formation his¬ 
tory of the host galaxy. Does this shift depend on the galaxy mass 
and does it therefore affect the relation between the GC metallicity 
and formation redshift? 

(viii) What is the reason for the existence of chemically enriched 
stellar populations in GCs? All present GC enrichment models do 
not rely on special, early-Universe physics, implying that the re¬ 
sponsible process should also take place in nearby YMCs. Because 
YMCs are seen to be devoid of gas after a few Myr, no evidence has 
yet been found for the multiple episodes of star formation required 
by the AGB and FRMS enrichment models. The search for ongoing 
or episodic star formation in YMCs should continue in order to test 
these scenarios. At the same time, an effort should be made to find 
evidence for the sweeping up of stellar ejecta by the protoplane¬ 
tary discs of pre-main sequence stars, which is a critical test for the 
EDA enrichment model. All of these tests can be made observation- 
ally at high spatial resolution in the local Universe. In addition, the 
evidence is mounting that none of the existing enrichment models 
can explain the observed abundance patterns (Bastian, Cabrera-Ziri 
& Salaris 2015), highlighting a clear need for new models, which 
(contrary to the previous generation) are not exclusively motivated 
by stellar evolution. 

(ix) Even if the majority of GCs has chemical abundance pat¬ 
terns consistent with a short formation time-scale, a handful of 
Galactic GCs show a spread or bimodality in their Ee abundance 
(e.g. ct;Cen). Some subset of these GCs likely did not form as reg¬ 
ular YMCs, but may be the remnants of dwarf galaxy nuclei (Lee 
et al. 2009). Given this example, which fraction of GCs formed 
through mechanisms other than regular YMC formation? Answer¬ 
ing this question requires current GC formation theories to be ex¬ 
tended such that the numbers of these ‘alternative’ GCs in present- 
day galaxy haloes can be predicted. Eor instance, coupling small- 
scale nuclear cluster formation models to the galaxy accretion his¬ 
tories of massive haloes from the Millennium Simulation enables 
their contribution to present-day GC populations to be quantified. A 
recent example of such an analysis was presented by Pfeffer et al. 
(2014), who find that at most a handful of Galactic GCs may be 
stripped dwarf galaxy nuclei. 

(x) In spite of all the evidence discussed in this work, it could 
be possible that for some reason the GCs observed today are not 
the products of regular cluster formation at z > 2. However, we 
have also seen that the products of regular YMC formation in the 
high-redshift Universe should have survived until the present day, 
and that they should have properties very similar to those of the GC 
populations seen in nearby galaxies. If GCs formed differently, the 
question arises what happened to large numbers of YMCs that must 
have formed across cosmic time - in the high-pressure environ¬ 
ments of z > 2 galaxies, a single YMC should have been formed 
for every ~ 10® M© in stars. While our two-phase model requires 
a large fraction of these to have been destroyed, it seems unlikely 
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that they could have disappeared altogether, as no known mecha¬ 
nism exists for disrupting a cluster with a mass of 10® M© and a 
radius of 1.5 pc within a Hubble time. Even in the [Fe/H] = —0.6 
galaxy model, which is the most disruptive environment consid¬ 
ered here, the short dynamical time of such a cluster allows it to 
respond adiabatically to tidal perturbations, leading to a disruption 
time-scale of ~ 10^ Gyr. Next to having to reproduce the observed 
properties of nearby GC populations, alternative GC formation the¬ 
ories would therefore also face the challenge of explaining why no 
products of regular YMC formation in z > 2 galaxies would have 
survived until the present day. 

The above list of open questions represents a relatively mod¬ 
est sample from the wide range of unsolved problems in the field 
of GC formation and evolution. Rather than being exhaustive, it 
is intended to highlight a number of areas that would particularly 
benefit from a concentrated effort in the next several years. As in¬ 
dicated above, many of these questions can be addressed using the 
new observational facilities that should come online in the near fu¬ 
ture (e.g. JWST) or have recently started operations (e.g. ALMA). 
In terms of theoretical or numerical work, several of the above 
questions can be addressed by considering specific physical mecha¬ 
nisms in isolation, whereas others require a more holistic approach. 
In order to test GC formation theories, it will be essential to de¬ 
velop a direct connection between GC formation physics and the 
observable properties of the GC population at z = 0. For instance, 
the two-phase GC formation model presented here is highly suit¬ 
able for implementation in semi-analytic models of galaxy forma¬ 
tion, as well as for verification and improvement using numerical 
simulations of galaxy formation. Above all, the current state of the 
field unequivocally shows that GC formation is a multi-scale, multi¬ 
physics problem, that is being constrained from a wide range of dif¬ 
ferent approaches, both in terms of the physics considered and the 
methods used. The only way of turning the physics of GC forma¬ 
tion into a quantitative and predictive science is to keep combining 
these different insights. 


7 SUMMARY AND CONCLUSIONS 

The purpose of this work is to address the question whether the 
products of young massive cluster (YMC, M > 10^ M©) for¬ 
mation in the high-redshift Universe may have survived until the 
present day, and if these relics would be consistent with the proper¬ 
ties of the globular cluster (GC) populations seen in nearby galax¬ 
ies. In §2, we summarised the state of the literature related to this 
question. Using the observational (§2.1) and theoretical (§2.2) con¬ 
straints on GC formation from nearby YMC formation, star for¬ 
mation in > 2 galaxies, and present-day GC populations, we 
derived a simple but quantitative, two-phase 'shaken, then stirred' 
model for GC formation (§3). Its main conclusions are as follows. 

(i) The properties of GCs are consistent with being the prod¬ 
ucts of the regular star and cluster formation process in the high- 
pressure environments of z > 2 galaxies, i.e. at the epoch of the 
peak cosmic SFR density. The early evolution of these clusters 
was characterised by rapid tidal disruption through encounters with 
dense gas clouds within the gas-rich host galaxy discs. However, 
the clusters were saved from tidal destruction in their dense natal 
environments by the redistribution of matter during galaxy merg¬ 
ers and/or tidal stripping, which allowed them to be ‘liberated’ into 
their host galaxy haloes, leading to their spatial distribution as ob¬ 


served today and marking the onset of an episode characterised by 
slow tidal disruption. 

(ii) This simple model reproduces a wide range of observables at 
z = 0 across a broad range of galaxies (§4): (1) the mass spectrum 
of GCs, (2) the relation between the number of GCs per unit stellar 
mass and the metallicity or galaxy stellar mass, (3) the colour and 
metallicity bimodality, and (4) the constant ratio between the GC 
population mass and the host galaxy’s dark matter halo mass. This 
variety of successes of the model is compelling, as it shows that 
combining the physics of star cluster formation and hierarchical 
galaxy formation naturally leads to the GC populations observed at 
z = 0. 

(iii) We find that including the rapid-disruption phase in GC for¬ 
mation models is essential to avoid making several ad-hoc assump¬ 
tions that would be necessary when only capturing the evolution of 
GCs in galaxy haloes. In halo-only models, the relations between 
the number of GCs and their host galaxies are attained after an en¬ 
tire Hubble time of GC evolution by evaporation. By contrast, the 
new model predicts that the relations between the number of GCs 
and their host galaxies were already in place at z = 1-2, whereas 
(primarily) the radii (and to a limited degree the masses) of GCs 
were shaped by their quiescent evolution in galaxy haloes. 

(iv) The high masses of GCs are a remnant of the Toomre mass 
during their formation, requiring surface densities of Y = 100- 
300 M© pc“^ and volume densities of pisM = 1-5 M© pc“®, 
which were common conditions in galaxies at z > 2. Due to the 
disruptive, high-pressure environments in these galaxies, only mas¬ 
sive clusters could survive until the time of the first galaxy merger 
and the corresponding migration of clusters into the galaxy halo 
(i.e. when the rapid-disruption phase ended and the clusters became 
GCs according to the definition of §1). For this reason, GC are typ¬ 
ically massive. 

(v) As the cosmic density decreases with cosmic time, so do 
both the Toomre mass and the galaxy merger rate. As a result, at 
lower redshifts fewer new-born clusters can survive long enough 
to escape into galaxy haloes, making it highly unlikely for them 
to survive over cosmic time-scales. For this reason, GCs are typi¬ 
cally old. Because the metallicity increases with cosmic time and 
galaxy mass, the peaking galaxy merger and GC survival rates at 
high redshift and low galaxy masses also imply that GCs are typi¬ 
cally metal-poor. 

(vi) During the hierarchical growth of galaxies, the native GC 
population remains associated with the host galaxy spheroid, 
whereas accreted GCs from stripped satellite galaxies mainly re¬ 
side in the halo. Because even the native GCs must have escaped 
the star-forming component of the host galaxy to survive over long 
time-scales, the typical galactocentric radii occupied by GCs ex¬ 
tend further than those of field stars, showing that a substantial 
fraction of the GC population resides in the halo. 

(vii) All models for the formation of an enriched population of 
stars in GCs predict that enrichment requires a minimum cluster 
mass (or density). By comparing the enriched stellar mass in GCs 
and in the field to the mass loss histories predicted by our model, 
we find that this minimum mass-scale is M r\j 10®-® M©. If true, 
the high initial masses of GCs (Mi ^ 10® M©, even for a current 
GC mass of M ~ 10^ M©) ensure that the minimum mass criterion 
is generally satisfied and hence the majority of GCs should display 
the observed multiple populations. 

In §1, we defined a GC as “a gravitationally-bound, stellar cluster 
that in terms of its position and velocity vectors does not coincide 
with the presently star-forming component of its host galaxy”. This 
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definition excludes clusters occupying the same part of position- 
velocity space as the host galaxy’s gaseous component. It thereby 
minimises GC destruction and thus represents a definition anal¬ 
ogous to an ‘anthropic principle’ for long-lived stellar clusters. 
We preferred this definition over previous ones adopted in the lit¬ 
erature, based on metallicity (‘metal-poor’), mass (M = 10^- 
10® M©), age (r ~ 10^® yr), location (‘in the halo’), or chemical 
abundance variations (‘multiple stellar populations’). The results 
summarised above show that all of these criteria are implied by 
combining our definition with the physics of GC formation. 

In the presented model, the ‘universal’ properties of GC pop¬ 
ulations expressed by the above range of criteria emerges irrespec¬ 
tively of their widely-varying present-day environments because 
the conditions of their formation were so similar. However, con¬ 
trary to previous GC formation theories, our model also allows ex¬ 
ceptions to these criteria, e.g. due to the stochastic nature of galaxy 
mergers, which in individual cases will occur on time-scales shorter 
or longer than the cosmic average. For instance, if a merger takes 
place immediately after a burst of massive cluster formation, the 
clusters are liberated into the halo shortly after their birth, giving 
rise to (possibly metal-rich) lower-mass, young GCs. 

The two-phase GC formation model explains a substantial 
number of observed characteristics of GCs (§4) and has a broad 
range of associated predictions and implications (§5), e.g. for con¬ 
straining the physical processes involved in GC formation (§5.1- 
5.12), as well as for testing the leading models for multiple stel¬ 
lar populations in GCs (§5.13-5.15). For instance, it is found that 
GCs were typically only 2-3 times more massive at birth and that 
most GCs could not have survived the high gas densities required 
for pristine gas accretion. These findings challenge the AGB and 
FRMS enrichment models. 

Next to these insights gained from the model, several open 
questions remain. These have been discussed in §6, but one par¬ 
ticular question is reiterated here. The evidence discussed in this 
work suggests that the products of regular YMC formation in the 
high-redshift Universe should have survived until the present day, 
and that they should have properties very similar to those of the GC 
populations seen in nearby galaxies. This conclusion lends support 
to the idea that YMCs in nearby galaxies can be used to test GC for¬ 
mation models. If GCs would not be the products of regular cluster 
formation at z > 2, the field would be facing a new problem: given 
the enormous number of YMCs produced in the history of the Uni¬ 
verse and their known longevity, where did they go? 

The ideal cosmic conditions for GC formation existed at z > 2 
due to the high mass-scales for gravitational instability and the fre¬ 
quent galaxy mergers that redistributed GCs into the silent waters 
of galaxy haloes. The likelihood of GC formation has since de¬ 
clined, but GCs can still form in the local Universe. When young, 
massive clusters are liberated into their host galaxy halo by a 
merger, their long-term survival has been secured, effectively turn¬ 
ing them into GCs. Although the GC population does erode with 
time, it presently does so at a much lower rate than in the past - the 
cluster disruption rate peaked in the host galaxy discs where GCs 
formed. Once thought to be the relics of a bygone era, the evidence 
discussed in this work suggests that GC populations will survive, 
and in rare cases even grow, for at least another Hubble time. 

ACKNOWLEDGEMENTS 

JMDK is greatly indebted to Nate Bastian, Jim Dale, Bruce 
Elmegreen, Mark Gieles, Henny Earners, Chervin Laporte, Spren 


Larsen, Steve Longmore, Simon Portegies Zwart, and Marina Rej- 
kuba for several years of discussions from which large parts of this 
work were drawn. In addition, Nate Bastian, Bill Harris, Gretchen 
Harris, and Steve Longmore are gratefully acknowledged for their 
careful reading of the manuscript. JMDK thanks Angela Adamo, 
Shy Genel, Michael Hilker, Phil Hopkins, Mark Krumholz, Clau¬ 
dia Lagos, Ben Moster, Thorsten Naab, Francois Schweizer, Linda 
Tacconi, and Simon White for helpful conversations and sugges¬ 
tions. The anonymous referee is thanked for constructive comments 
that improved the presentation of this work. 


REFERENCES 

Adamo A., Kruijssen J. M. D., Bastian N., Silva-Villa E., Ryon J., 2015, 
MNRAS, 452, 246 

Adams F. C., Hollenbach D., Laughlin G., Gorti U., 2004, ApJ, 611, 360 
Anders P., de Grijs R., Fritze-v. Alvensleben U., Bissantz N., 2004, 
MNRAS, 347, 17 

Ashman K. M., Zepf S. E., 1992, ApJ, 384, 50 
Bastian N., 2008, MNRAS, 390, 759 
Bastian N. et al., 2012, MNRAS, 419, 2606 

Bastian N., Cabrera-Ziri L, Davies B., Larsen S. S., 2013a, MNRAS, 436, 
2852 

Bastian N., Cabrera-Ziri L, Salaris M., 2015, MNRAS, 449, 3333 
Bastian N., Ercolano B., Gieles M., Rosolowsky E., Scheepmaker R. A., 
Gutermuth R., Efremov Y, 2007, MNRAS, 379, 1302 
Bastian N., Earners H. J. G. L. M., de Mink S. E., Longmore S. N., Good¬ 
win S. R, Gieles M., 2013b, MNRAS, 436, 2398 
Bastian N., Saglia R. R, Goudfrooij R, Kissler-Patig M., Maraston C., 
Schweizer K, Zoccali M., 2006, A&A, 448, 881 
Bastian N., Tranche G., Konstantopoulos 1. S., Miller B. W., 2009, ApJ, 
701,607 

Baumgardt H., 2001, MNRAS, 325, 1323 
Baumgardt H., Makino J., 2003, MNRAS, 340, 227 
Baumgardt H., Parmentier G., Gieles M., Vesperini E., 2010, MNRAS, 
401, 1832 

Bekki K., 2006, MNRAS, 367, L24 

Bell C. P. M., Naylor T, Mayne N. J., Jeffries R. D., Littlefair S. R, 2013, 
MNRAS, 434, 806 
Bell E. F. et ah, 2008, ApJ, 680, 295 
Bigiel F. et ah, 2011, ApJ, 730, L13 

Boley A. C., Lake G., Read J., Teyssier R., 2009, ApJ, 706, L192 

Bondi H., 1952, MNRAS, 112, 195 

Bondi H., Hoyle E, 1944, MNRAS, 104, 273 

Bonnell 1. A., Clark R, Bate M. R., 2008, MNRAS, 389, 1556 

Bouche N. et ah, 2010, ApJ, 718, 1001 

Bournaud R, Elmegreen B. G., Elmegreen D. M., 2007, ApJ, 670, 237 

Brodie J. P. et ah, 2014, ApJ, 796, 52 

Brodie J. R, Strader J., 2006, ARA&A, 44, 193 

Burkert A., Tremaine S., 2010, ApJ, 720, 516 

Cabrera-Ziri L, Bastian N., Davies B., Magris G., Bruzual G., Schweizer 
R, 2014, MNRAS, 441, 2754 

Carraro G., Villanova S., Demarque R, McSwain M. V, Piotto G., Bedin 
L. R.,2006, ApJ, 643, 1151 

Carretta E., Bragaglia A., Gratton R. G., Recio-Blanco A., Lucatello S., 
D’Orazi V., Cassisi S., 2010, A&A, 516, A55 
Chernoff D. R, Weinberg M. D., 1990, ApJ, 351, 121 
Clarkson W. L, Ghez A. M., Morris M. R., Lu J. R., Stolte A., McCrady 
N., Do T, Yelda S., 2012, ApJ, 751, 132 
Conroy C., 2012, ApJ, 758, 21 
Conroy C., Loeb A., Spergel D. N., 2011, ApJ, 741, 72 
Conroy C., Spergel D. N., 2011, ApJ, 726, 36 
Cook D. O. et ah, 2012, ApJ, 751, 100 

Cottaar M., Meyer M. R., Andersen M., Espinoza R, 2012, A&A, 539, 
A5 

Cowie L. L., Songaila A., Hu E. M., Cohen J. G., 1996, AJ, 112, 839 


© 2015 RAS, MNRAS 000 , 1-29 


28 J. M. D. Kruijssen 


Daddi E. et aL, 2010, ApJ, 714, LI 18 
Danielson A. L. R. et al, 2013, MNRAS, 436, 2793 
de Juan Ovelar M., Kruijssen J. M. D., Bressert E., Testi L., Bastian N., 
Canovas H., 2012, A&A, 546, LI 

de Mink S. E., Pols O. R., Langer N., Izzard R. G., 2009, A&A, 507, LI 
Decressin T., Baumgardt H., Kroupa R, 2008, A&A, 492, 101 
Decressin T., Meynet G., Charbonnel C., Prantzos N., Ekstrom S., 2007, 
A&A, 464, 1029 

D’Ercole A., Vesperini E., D’Antona R, McMillan S. L. W., Recchi S., 
2008, MNRAS, 1228 

Dinescu D. L, Girard T. M., van Altena W. K, 1999, AJ, 117, 1792 
Downes D., Solomon R M., 1998, ApJ, 507, 615 
Durrell R R. et al, 2014, ApJ, 794, 103 
Elmegreen B. G., 1983, MNRAS, 203, 1011 
Elmegreen B. G., 2000, ApJ, 530, 277 
Elmegreen B. G., 2008, ApJ, 672, 1006 
Elmegreen B. G., 2010, ApJ, 712, L184 
Elmegreen B. G., Efremov Y. N., 1997, ApJ, 480, 235 
Elmegreen B. G., Elmegreen D. M., 2001, AJ, 121, 1507 
Elmegreen B. G., Elmegreen D. M., 2005, ApJ, 627, 632 
Elmegreen B. G., Falgarone E., 1996, ApJ, 471, 816 
Elmegreen B. G., Hunter D. A., 2010, ApJ, 712, 604 
Erb D. K., Steidel C. C., Shapley A. E., Pettini M., Reddy N. A., Adel- 
berger K. L., 2006, ApJ, 646, 107 
Evans, IIN. J. et al, 2009, ApJS, 181, 321 

Fakhouri O., Ma C.-R, Boylan-Kolchin M., 2010, MNRAS, 406, 2267 

Fall S. M., Rees M. J., 1977, MNRAS, 181, 37P 

Fall S. M., Rees M. J., 1985, ApJ, 298, 18 

Fall S. M., Zhang Q., 2001, ApJ, 561, 751 

Ferrarese L., 2002, ApJ, 578, 90 

Ferrarese L. et al., 2012, ApJS, 200, 4 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Eiger D. R, McLean 1. S., Morris M., 1999, ApJ, 514, 202 

Eiger D. R et al, 2002, ApJ, 581, 258 

Fisher D. B. et al, 2014, ApJ, 790, L30 

Forbes D. A., Bridges T., 2010, MNRAS, 404, 1203 

Forbes D. A., Forte J. C., 2001, MNRAS, 322, 257 

Forster Schreiber N. M. et ah, 2009, ApJ, 706, 1364 

Frenk C. S., White S. D. M., 1980, MNRAS, 193, 295 

Fuji! M. S., Saitoh T. R., Portegies Zwart S. R, 2012, ApJ, 753, 85 

Fukushige T., Heggie D. C., 2000, MNRAS, 318, 753 

Gebhardt K. et al, 2000, ApJ, 539, L13 

Genel S., Genzel R., Bouche N., Naab T., Sternberg A., 2009, ApJ, 701, 
2002 

Genzel R. et ah, 2011, ApJ, 733, 101 
Genzel R. et ah, 2010, MNRAS, 407, 2091 

Georgiev 1. Y, Puzia T. H., Goudfrooij R, Hilker M., 2010, MNRAS, 
406, 1967 

Gieles M., 2009, MNRAS, 394, 2113 

Gieles M., Bastian N., Earners H. J. G. L. M., Mout J. N., 2005, A&A, 
441, 949 

Gieles M., Baumgardt H., 2008, MNRAS, 389, L28 
Gieles M., Heggie D. C., Zhao H., 2011, MNRAS, 413, 2509 
Gieles M., Larsen S. S., Scheepmaker R. A., Bastian N., Haas M. R., 
Earners H. J. G. L. M., 2006a, A&A, 446, L9 
Gieles M., Portegies Zwart S. R, Baumgardt H., Athanassoula E., Earners 
H. J. G. L. M., Sipior M., Leenaarts J., 2006b, MNRAS, 371, 793 
Girichidis R, Federrath C., Allison R., Banerjee R., Klessen R. S., 2012, 
MNRAS, 420, 3264 

Gnedin O. Y, Hernquist L., Ostriker J. R, 1999, ApJ, 514, 109 

Gnedin O. Y, Ostriker J. R, 1997, ApJ, 474, 223 

Goddard Q. E., Bastian N., Kennicutt R. C., 2010, MNRAS, 405, 857 

Gouliermis D. A., Hony S., Klessen R. S., 2014, MNRAS, 439, 3775 

Gratton R. G., Carretta E., Bragaglia A., 2012, A&A Rev., 20, 50 

Haisch, Jr. K. E., Lada E. A., Lada C. J., 2001, ApJ, 553, L153 

Harris G. L. H., Harris W. E., 2011, MNRAS, 410, 2347 

Harris W. E., 1991, ARA&A, 29, 543 

Harris W. E., 1996, AJ, 112, 1487 


Harris W. E., Harris G. L., Hudson M. J., 2015, ApJ, 806, 36 
Harris W. E., Harris G. L. H., 2002, AJ, 123, 3108 
Harris W. E., Harris G. L. H., Alessi M., 2013, ApJ, 772, 82 
Hartmann M., Debattista V. R, Seth A., Cappellari M., Quinn T. R., 2011, 
MNRAS, 418, 2697 

Henault-Brunet V. et al., 2012, A&A, 546, A73 

Hodge J. A., Carilli C. L., Walter R, de Blok W. J. G., Riechers D., Daddi 
E., Lentati L., 2012, ApJ, 760, 11 
Holtzman J. A. et al., 1992, AJ, 103, 691 
Hopkins A. M., Beacom J. R, 2006, ApJ, 651, 142 
Hopkins P. R, Cox T. J., Younger J. D., Hernquist L., 2009, ApJ, 691, 
1168 

Hoyle R, Lyttleton R. A., 1941, MNRAS, 101, 227 
Hudson M. J., Harris G. L., Harris W. E., 2014, ApJ, 787, L5 
Jordan A. et al., 2007, ApJS, 171, 101 
Katz H., Ricotti M., 2013, MNRAS, 432, 3250 
Katz H., Ricotti M., 2014, MNRAS, 444, 2377 
Kennicutt, Jr. R. C., 1998, ApJ, 498, 541 
Kim W.-T., Ostriker E. C., 2001, ApJ, 559, 70 
Kormendy J., Ho L. C., 2013, ARA&A, 51, 511 
Krause M., Charbonnel C., Decressin T., Meynet G., Prantzos N., 2013, 
A&A, 552, A121 

Kravtsov A. V., Gnedin O. Y, 2005, ApJ, 623, 650 
Kruijssen J. M. D., 2011, PhD Thesis, Utrecht University, Formation and 
evolution of star clusters and their host galaxies 
Kruijssen J. M. D., 2012, MNRAS, 426, 3008 
Kruijssen J. M. D., 2013, arXiv: 1304.4600 

Kruijssen J. M. D., 2014, Classical and Quantum Gravity, 31, 244006 
Kruijssen J. M. D., Cooper A. R, 2012, MNRAS, 420, 340 
Kruijssen J. M. D., Longmore S. N., 2013, MNRAS, 435, 2598 
Kruijssen J. M. D., Longmore S. N., Elmegreen B. G., Murray N., Bally 
J., Testi L., Kennicutt R. C., 2014, MNRAS, 440, 3370 
Kruijssen J. M. D., Maschberger T., Moeckel N., Clarke C. J., Bastian N., 
Bonnell 1. A., 2012a, MNRAS, 419, 841 
Kruijssen J. M. D., Mieske S., 2009, A&A, 500, 785 
Kruijssen J. M. D., Pelupessy F. I., Earners H. J. G. L. M., Portegies Zwart 
S. R, Bastian N., Icke V., 2012b, MNRAS, 421, 1927 
Kruijssen J. M. D., Pelupessy F. I., Earners H. J. G. L. M., Portegies Zwart 
S. R, Icke V., 2011, MNRAS, 414, 1339 
Kruijssen J. M. D., Portegies Zwart S. R, 2009, ApJ, 698, L158 
Krumholz M. R., Klein R. I., McKee C. R, 2012, ApJ, 754, 71 
Krumholz M. R., McKee C. R, 2005, ApJ, 630, 250 
Lada C. J., Lada E. A., 2003, ARA&A, 41, 57 

Earners H. J. G. L. M., Baumgardt H., Gieles M., 2010, MNRAS, 409, 
305 

Earners H. J. G. L. M., Gieles M., 2006, A&A, 455, L17 
Earners H. J. G. L. M., Gieles M., Bastian N., Baumgardt H., Kharchenko 
N. V., Portegies Zwart S., 2005, A&A, 441, 117 
Earners H. J. G. L. M., Gieles M., Portegies Zwart S. R, 2005, A&A, 429, 
173 

Earners H. J. G. L. M., Kruijssen J. M. D., Bastian N., Rejkuba M., Hilker 
M., Kissler-Patig M., 2015, in preparation 
Larsen S. S., 2000, MNRAS, 319, 893 
Larsen S. S., 2004, A&A, 416, 537 
Larsen S. S., 2009, A&A, 494, 539 

Larsen S. S., Brodie J. R, Forbes D. A., Strader J., 2014, A&A, 565, A98 

Larsen S. S., Strader J., Brodie J. R, 2012, A&A, 544, L14 

Lee J.-W., Kang Y.-W., Lee J., Lee Y.-W., 2009, Nature, 462, 480 

Li H., Gnedin O. Y, 2014, ApJ, 796, 10 

Longmore S. N. et al., 2014, Protostars and Planets VI, 291 

Mackey A. D., Gilmore G. R, 2003, MNRAS, 340, 175 

Mackey A. D., van den Bergh S., 2005, MNRAS, 360, 631 

Madrid J. R, Hurley J. R., Sippel A. C., 2012, ApJ, 756, 167 

Magorrian J. et al., 1998, AJ, 115, 2285 

Mannucci F. et al., 2009, MNRAS, 398, 1915 

Maschberger T., Clarke C. J., Bonnell I. A., Kroupa R, 2010, MNRAS, 
404, 1061 


© 2015 RAS, MNRAS 000 , 1-29 


The origin of globular clusters 29 


Maxwell A. J., Wadsley J., Couchman H. M. R, Sills A., 2014, MNRAS, 
439, 2043 

McLaughlin D. E., Fall S. M., 2008, ApJ, 679, 1272 
Moore B., Diemand J., Madau R, Zemp M., Stadel J., 2006, MNRAS, 
368, 563 

Moster B. R, Naab T., White S. D. M., 2013, MNRAS, 428, 3121 
Muratov A. L., Gnedin O. Y, 2010, ApJ, 718, 1266 
Offner S. S. R., Klein R. L, McKee C. F., Krumholz M. R., 2009, ApJ, 
703, 131 

Ostriker J. R, Spitzer L. J., Chevalier R. A., 1972, ApJ, 176, L51+ 

Rarmentier G., Gilmore G., 2007, MNRAS, 377, 352 

Reebles R J. E., 1984, ApJ, 277, 470 

Reebles R J. E., Dicke R. H., 1968, ApJ, 154, 891 

Reng E. W et ah, 2006, ApJ, 639, 95 

Reng E. W. et ah, 2008, ApJ, 681, 197 

Reng Y, Maiolino R., Cochrane R., 2015, Nature, 521, 192 

Reng Y-j. et ah, 2010, ApJ, 721, 193 

Rfeffer J., Griffen B. E, Baumgardt H., Milker M., 2014, MNRAS, 444, 
3670 

Rlanck Collaboration et aL, 2014, A&A, 571, A16 
Rontzen A., Governato R, 2012, MNRAS, 421, 3464 
Rortegies Zwart S. R, Hut R, Makino J., McMillan S. L. W., 1998, A&A, 
337, 363 

Rortegies Zwart S. R, McMillan S. L. W, Gieles M., 2010, ARA&A, 48, 
431 

Rrieto J. L., Gnedin O. Y, 2008, ApJ, 689, 919 
Rathborne J. M. et aL, 2015, ApJ, 802, 125 
Renaud R, Gieles M., 2013, MNRAS, 431, L83 

Rieder S., Ishiyama T., Langelaan R, Makino J., McMillan S. L. W., Rorte¬ 
gies Zwart S., 2013, MNRAS, 436, 3695 
Rochau B., Brandner W., Stolte A., Gennaro M., Gouliermis D., Da Rio 
N., Dzyurkevich N., Henning T., 2010, ApJ, 716, L90 
Rodighiero G. et aL, 2011, ApJ, 739, L40 

Rosotti G. R, Dale J. E., de Juan Ovelar M., Rubber D. A., Kruijssen 
J. M. D., Ercolano B., Walch S., 2014, MNRAS, 441, 2094 
Sabbi E. et aL, 2012, ApJ, 754, L37 
Schechter R, 1976, ApJ, 203, 297 
Schruba A. et aL, 2011, AJ, 142, 37 
Schweizer R, 1982, ApJ, 252, 455 

Schweizer R, 1987, in Nearly Normal Galaxies. From the Rlanck Time to 
the Rresent, New York, Springer-Verlag, S. M. Faber, ed., pp. 18-25 
Schweizer R, Miller B. W., Whitmore B. C., Fall S. M., 1996, AJ, 112, 
1839 

Schweizer R, Seitzer R, 1998, AJ, 116, 2206 
Searle L., Zinn R., 1978, ApJ, 225, 357 

Shapiro K. L., Genzel R., Forster Schreiber N. M., 2010, MNRAS, 403, 
L36 

Sijacki D., Vogelsberger M., Genel S., Springel V., Torrey R, Snyder G. R, 
Nelson D., Hernquist L., 2015, MNRAS, 452, 575 
Silva-Villa E., Adamo A., Bastian N., 2013, MNRAS, 436, L69 
Spider L. R., Forbes D. A., 2009, MNRAS, 392, LI 
Spitzer L., 1987, Dynamical evolution of globular clusters. Rrinceton, NJ, 
Rrinceton University Rress, 1987, 191 p. 

Spitzer, Jr. L., 1958, ApJ, 127, 17 
Springel V. et aL, 2005, Nature, 435, 629 
Swinbank A. M. et aL, 2011, ApJ, 742, 11 

Swinbank A. M., Small L, Sobral D., Theuns T., Best R. N., Geach J. E., 
2012, ApJ, 760, 130 

Tacconi L. J. et aL, 2010, Nature, 463, 781 

Tacconi L. J. et aL, 2013, ApJ, 768, 74 

Tilley D. A., Rudritz R. E., 2007, MNRAS, 382, 73 

Tonini C., 2013, ApJ, 762, 39 

Toomre A., 1964, ApJ, 139, 1217 

van den Bosch F. C., Aquino D., Yang X., Mo H. J., Rasquali A., McIntosh 
D. H., Weinmann S. M., Kang X., 2008, MNRAS, 387, 79 
VandenBerg D. A., Brogaard K., Leaman R., Casagrande L., 2013, ApJ, 
775,134 

Vesperini E., 2001, MNRAS, 322, 247 


Vesperini E., Heggie D. C., 1997, MNRAS, 289, 898 
Vesperini E., McMillan S. L. W., D’Antona R, D’Ercole A., 2010, ApJ, 
718,L112 

Vesperini E., Zepf S. E., Kundu A., Ashman K. M., 2003, ApJ, 593, 760 
Walker A. R. et aL, 2011, MNRAS, 415, 643 

Walker D., Longmore S. N., Bastian N., Kruijssen J. M. D., Rathborne 
J. M., Jackson J. M., Foster J. B., Contreras Y, 2015, MNRAS, 449, 
715 

Walker M. G., Renarrubia J., 2011, ApJ, 742, 20 
Weinberg M. D., 1994, AJ, 108, 1403 
White S. D. M., Frenk C. S., 1991, ApJ, 379, 52 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
Whitmore B. C., Zhang Q., Leitherer C., Fall S. M., Schweizer R, Miller 
B. W., 1999, AJ, 118, 1551 
Wright E. L., 2006, RASP, 118, 1711 

Wright N. J., Barker R. J., Goodwin S. R, Drake J. J., 2014, MNRAS, 
438, 639 

Zhang Q., Fall S. M., 1999, ApJ, 527, L81 
Zinn R., 1985, ApJ, 293, 424 

This paper has been typeset from a TgX/ UT^ file prepared by the 
author. 


© 2015 RAS, MNRAS 000 , 1-29 


